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I.  Introduction 


Spinning  projectiles  containing  liquid  payloads  have  shown  spectacular  instabilities  in 
flight.  The  first  theoretical  work  by  Stewartson^  in  1959  showed  that  a  spinning  inviscid 
liquid  could  exert  a  very  large  side  moment  for  certain  values  of  r,  the  ratio  of  coning  rate 
to  spin  rate  (symbols  are  defined  in  the  List  of  Symbols  near  the  back  of  this  report).  In 
19G5  Wedemeyer^  added  a  boundary  layer  on  the  container  walls;  his  predicted  values  of 
the  resonant  side  moment  were  verified  by  experiments  using  gyroscopes.  This  Stewartson- 
Wodemejer  (SW)  boundary- layer  theory  was  originally  limited  to  pressure  moments,  fully- 
or  partially-filled  cylindrical  cavities,  and  fully  spun-up  liquids.  Later  authors!3-*1  extended 
the  theory  to  include  centred  rods,  viscous  moments,  partially  spun-up  liquids,  two  immis¬ 
cible  liquids,  and  fully- filled  spheroidal  cavities.  In  addition,  a  relation  between  liquid 
side  moment  and  liquid  roll  moment  was  derived!91  to  justify  experimental  observation  of 
a  correlation  between  these  moments. 

In  197S  Kitchens,  Gerber  and  Sedney!101  developed  a  hybrid  linear  Navier-Stokes  the¬ 
ory  (hereafter  called  KGS)  that  satisfied  the  boundary  conditions  on  the  cylindrical  wall 
exactly  by  solving  a  sixth-order  system  of  ordinary  differential  equations,  but  it  required 
the  addition  of  boundary  layers  on  the  end-walls  to  satisfy  the  no-slip  condition  on  these 
walls.  This  hybrid  high-Reynolds-number  theory  has  been  extended  to  predict  the  com¬ 
plete  moment  exerted  by  both  fully  spun-up  liquids111’*21  and  partially  spun-up  liquids.!13,14! 

For  very  low  Reynolds  numbers  (Re  <  100),  Vaughn  et  al.!1*1,  Strikwerda  et  al. 118,171 
and  Rosenblat  et  al.11*1  have  developed  Computational  Fluid  Dynamics  (CFD)  codes  to 
predict  moments  for  a  fully-filled  cylinder  coning  at  a  constant  angle.  These  predic¬ 
tions  showed  good  agreement  with  experiments,  but  this  CFD  approach  requires  consider¬ 
able  time  on  large  computing  machines.  Herbert!19,201  developed  an  approximate  spectral 
method  that  requires  rather  modest  computational  effort  and  gives  good  agreement  with 
the  more  precise  CFD  results  for  Rt  <  100. 

Hall,  Sedney  and  Gerber*211  have  recently  modified  the  hybrid  Navier-Stokes  theory 
to  eliminate  the  approximation  of  the  end-wall  boundary  conditions.  This  was  done  by  re¬ 
placing  the  simple  Stewartson-Wedemeyer  spatial  eigenvalues  by  special  eigenvalues  whose 
eigenfunctions  can  be  combined  to  satisfy  all  the  boundary  conditions.  Once  a  table 
of  eigenvalues  has  been  constructed,  this  HSG  method  can  compute  liquid  moments  for 
Re  <  2500  on  VAX  8600-size  computers  in  less  than  five  minutes.  These  computations 
show  exce  .ent  agreement  with  the  CFD  calculations!221  for  Re  <  100  and  good  agreement 
with  the  boundary-layer  theories  for  Re  >  2000. 

The  KGS  and  HSB  theories  are  limited  to  a  fully-filled  cylindrical  container  performing 
constant-amplitude  coning  motion.  In  this  report,  we  will  extend  these  theories  to  partially- 
filled  cylinders  and  to  cylinders  with  a  central  rod.  We  will  also  allow  the  amplitude  of  the 
motion  to  change  slowly.  Next,  we  will  extend  the  theories  to  the  esse  of  two  immiscible 
liquids. 

In  the  implementation  of  the  SW,  KGS  and  HSG  theories,  the  liquid  moment  has 
usually  been  computed  by  integrating  the  pressure  and  wail  shears  over  the  surface  of  the 
cylindrical  container.  A  mathematically  equivalent  approach  is  to  integrate  the  angular 
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momentum  over  the  liquid  volume  and  then  differentiate  this  integral.  It  will  be  shown  that 
for  the  essentially  exact  HSG  theory,  the  two  methods  -  surface  and  volume  integration 
-  yield  the  same  value  for  the  liquid  moment.  However,  for  the  less  exact  3W  and  KGS 
theories,  the  volume  integral  method  gives  much  ‘better’  moment  values,  in  the  sense  that 
they  are  in  closer  agreement  with  HSG  results.  Indeed,  using  the  volume  integral,  the 
5W  theory  adequately  approximates  HSG  results  down  to  Re  =  500  and  KGS  adequately 
approximates  HSG  down  to  Re  =  1. 


II.  Equations  of  Liquid  Motion 


We  will  consider  a  projectile  with  a  cylindrical  cavity  of  radius  a  and  height  2c  and 
hence  of  fineness  ratio  A  =  c/a.  The  axis  of  the  cylinder  is  collinear  with  the  projectile’s 
axis  and  the  cylinder  center  of  mass  is  located  a  distance  h  from  the  projectile’s  center  of 
mass.  For  simplicity,  we  will  assume  that  h  =  0;  the  effect  of  a  nonzero  h  will  be  considered 
in  Appendix  A. 

R"  the  cavity  is  only  partially  filled,  if  the  liquid  is  fully  spun-up  and  if  the  centrifugal 
force  is  large  compared  to  the  aerodynamic  forces,  then  in  the  absence  of  coning  motion 
the  liquid  will  fill  the  space  between  the  outer  cylindrical  wall  and  an  inner  cylindrical  free 
surface  of  radius  b.  The  ratio  of  the  volume  of  this  inner  cylinder  to  the  volume  of  the 
payload  cavity  i3  b2/a2.  The  fill  ratio  for  the  payload  cavity  is  thus  1  —  b2 /a1. 

In  addi'/on  to  fully  and  partially  filled  cylindrical  payload  cavities,  the  case  of  a 
cylindrical  cavity  with  a  central  rod  of  radius  d  can  be  easily  treated  by  the  linear  theory. 
The  annular  region  between  the  inner  and  outer  cylindrical  walls  will  be  considered  fully 
filled  with  liquid. 

Although  the  differential  equations  of  motion  are  most  conveniently  expressed  in  earth- 
fixed  coordinates,  the  boundary  conditions  are  best  expressed  in  cylindrical  coordinates 
aligned  with  the  projectile’s  axis  of  symmetry.  Symmetry-axis  coordinates  have  one  free 
parameter:  angular  velocity  about  the  symmetry  axis.  If  this  angular  velocity  is  chosen 
to  be  the  projectile’s  spin  rate,  the  system  is  the  usual  missile-fixed  coordinate  system.  If 
this  angular  velocity  is  chosen  to  be  zero,  we  have  the  nonspinning  aeroballistic  system 
used  in  the  flight  mechanics  of  symmetric  missiles^. 

Earth-fixed  axes  .Ye,  Ye,  Zt  are  selected  so  that  X,  is  initially  along  the  velocity  vector1 
and  the  A'-axis  of  the  missile-aligned  system  will  be  along  the  projectile’s  axis  of  symmetry. 
The  angle  between  the  AT  and  AT,  axes  is  the  total  angle  of  attack,  at.  The  angle  between 
the  Zf- axis  and  the  normal  to  the  plane  formed  by  X  and  .Ye  will  be  called  <f>a,  where  4>a  is 
the  precessional  rate  of  the  angular  motion  performed  by  the  projectile  about  *he  velocity 
vector. 

Although  the  angular  motion  of  a  projectile  does  not  in  general  have  a  constant  pre¬ 
cessional  rate,  it  cam  be  shown  that  the  linear  motion  is  the  sum  of  two  circular  motions 
with  constant  precessional  rates.  The  objective  of  all  liquid  payload  theories  is  to  predict 

1  In  thia  report,  wm  mwim  that  the  velocity  vector  maintains  a  coact  ant  direction.  Th«  dbett  of  aerodynamic  forcaa  and 
gravity  arc  given  in  Reference  (23). 
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the  liquid  moment  response  to  constant-precession  angular  motion.  For  these  theories, 
projectile- aligned  coordinates  that  rotate  at  the  precessional  rate  (coning  aeroballistic  co¬ 
ordinates)  are  especially  convenient  since  the  boundary  conditions  are  then  independent 
of  time. 

Let  (ex,ev,ex)  be  unit  vectors  in  the  earth-fixed  system  and  let  (etc>e*e*«*c)  be  unit 
vectors  in  ihe  toning  aeroballistic  system.  These  unit  vectors  are  related  as  follows: 

exc  =  7«r  —  A'(ey  cos  <f>a  +  e,  sin  <j>a)  (1) 

eyc  =  Aet  +  7(evcos<k,  +  etsin<£a)  (2) 

e«  =  -ev3in^>a -(- ex  cos^o  (3) 


where 

K  —  sin  Oft,  7  *  cos  art 

The  angular  velocity  ft  of  the  coning  aeroballistic  coordinate  system  is  defined  by  the 
relation 

(erci  Cyc,  Cjc)  =  ft  x  (cre,  eye,  exc)  (4) 

It  can  be  easily  verified  that 

ft  =  t4>[£x  -  (<AV7)e,-]  (5) 

c*  *  7«*c  +  A'cye  (6) 

where 


•  <f>  is  the  projectile  spin  rate  relative  to  sin  inertial  frame; 

•  t  =  4>a/  <j>,  the  ratio  of  the  precessional  rate  to  the  spin  rate;  and 

•  K  —  KottT  .  Hence  =  K/K  and  we  see  that  e  is  a  measure  of  the  logarithmic 
damping  during  a  precession  cycle.  Zero  e  denotes  constant-amplitude  coning  motion. 

Cylindrical  coordinates  will  be  denoted  by  ( x,r,9 )  in  the  earth-fixed  system  and  by 
(x,  r,  rjj)  in  the  coning  aeroballistic  system.  (For  convenience,  we  will  assume  that  all  length 
variables  have  been  made  dimensionless  by  division  by  the  radius  a.)  From  Eqs.(l-3),  we 
have 


x  =  7X  f  A'  cos  %i>  (7) 

rcosd  =  f[cos(^  +  <t>a)  +  (7-  i) cos  ^cos^0]  -  xA' cos  <f>a  (8) 

rsintf  =  r[sin(V>  +  4>a)  +  (7  —  1)  cos  0  sin  da]  —  xA"  sin^tt  (9) 

For  small  K,  Eqs.(7-9)  reduce  to 

x  =*  x  +  ?R{Ke~'*}  (10) 

r  =  r  -  xfl{A'e“*}  (11) 

sin (0  -  VO  ®  -(x/r)R{iKe~'*}  (12) 
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where 


0  =  9  -  <t>a  =  9  —  r<f> 

<j>  =  j>t 
<t>a  =  4>at 

and  where  R{  }  denotes  the  real  part  of  a  complex  quantity. 

The  coning  cylindrical  components  of  the  velocity  of  a  point  on  the  projectile  are 
given  by 

i  =  0  (13) 

f  =  0  (14) 

=  r<j>(  1  -  r)  (15) 

The  earth-fixed  cylindrical  components  of  the  velocity  can  be  obtained  by  differentiating 
Eqs.(10-12): 

V,  =  a$R{i(f  —  l)r/C*}  (16) 

Vr  *  —aj>R{i(f  —  l)xA'*}  (17) 

V,  »  a^(r-i2{(/-l)xA:-}]  (18) 

where 

A'*  =*  AV*  *  AV*'*"** 

/  «  (l-u)r 

We  will  now  make  the  very  restrictive  assumption  that  the  liquid  is  in  steady-state 
response  to  the  coning  and  spinning  motion  of  the  projectile.  Theoretical  studies!3®!  have 
been  made  and  are  in  progress  to  determine  the  effect  of  partially  spun-up  liquid,  and  an 
experimental  study!24!  has  been  ciade  of  the  transient  response  to  coning  motion.  These 
studies  show  that  spin-up  and  cone-up  effects  are  large  and  important  to  a  complete 
understanding  of  the  liquid- payload  stability  problem. 

Nevertheless,  we  will  assume  that  the  liquid  velocity  components  and  liquid  pressure 
have  the  same  dependency  on  time  and  9  as  do  the  velocity  components  of  points  on  the 
projectile.  Accordingly,  we  introduce  three  dimensionless  complex  perturbation  velocities: 
ULtlii  lit  each  a  function  of  r  and  x,  and  write: 


V,  =»  -aiR{uiI<m) 

(19) 

Vr  =  -a^A{uA’*} 

(20) 

V,  -  aj>  (r  —  A{i>A'*}] 

(21) 

The  equilibrium  pressure  for  a  spinning  non-coning  cylinder  is 

•  fr2  -  r3 

„  0  (r0  <  r  <  1) 

(22) 

where 
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•  r0  is  either  b/a  or  d/a , 

•  pi  is  the  liquid  density. 


When  the  cylinder  is  performing  coning  motion,  the  outer  cylindrical  wall  is  located  at 
f  =  1.  The  equivalent  earth- fixed  radial  coordinate  of  the  wall  is 

rc  =  1  —  R{xK’}  (23) 

The  expression  for  the  earth-fixed  radial  coordinate  of  a  central  rod  is  similar  in  form  to 

(23): 

rr  =  r0-R{xKm)  (24) 

For  a  coning,  partially  filled  cylinder,  the  free  surface  of  the  liquid  is  a  slightly  perturbed 
version  of  a  cylindrical  surface  of  radius  r<j: 

ry-^l-S{2i(*)JC}]  (25) 

To  determine  note  that  the  free  surface  moves  with  the  radial  velocity  of  the  liquid  at 
r  =  r0.  Hence  differentiation  of  Eq.(25)  yields 


Thus 


Vr  =  -rnR{[i%(x)(af+  ~  V,/r0)  - 
„  /rx  »ii(ro,x) 

3»(I)  -  (TT7h 


The  complex  liquid  pressure  perturbation  function  £  is,  therefore,  defined  by  the  relation 
P-PD  +  PtaV  -/?{£*•}  ,  (rT,r,<r<rc)  (28) 

The  linearized  Navier-Stokes  equations  for  these  perturbation  functions  are 


d(ru)  ,  ota 

—  -,u+r  a7 


■it  ,  1  ,  ifo  2(a - >a) 

*(/-l)u-2u  *  “ar  +  J?e  [dr3  +  dx3  +  r  dr  ~  r3 

,  «  *E  .  1  ,  I5ji  2i(u  —  *ii) 

■(/  -  i)k+2u  =  f +  sla^  +  5?  +  ;S - ? — 

*  _  0£  1  [d3^  d3^  1  dia  a:] 

*(/  l)ut  -  dx  +  Re  dr3  +  dx3  +  r  dr  ~  r3 


Each  of  the  four  perturbation  functions  can  now  be  assumed  to  be  the  sum  of  products  of 
two  variables:  the  first  a  function  of  r  alone  and  the  second  a  function  of  x  alone: 

la  «  *(1  -  f)rHo(x)  -  £  Wk(r)Hk(x)  (33) 

kmO 

u  =  -  hlTfl  G„(i)  +  E .i.(r)C,(z)  (34) 
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(35) 


Go(x)  +  Y,  Vk(r)Gk(x) 
i=0 

p  =  -(1  -  f)2>'Go(x)  +  YPk(r)Gk(x)  (36) 

k= 0 

In  References  [4,  10-12], 

for  A0  =  0,  Ho  —  1  and  Go  =  x  (37) 

for  A*  7^  0,  Hk  —  cos(Akx)  and  Gi,  =  sin(Akx)  (38) 

while  in  HSG,  Reference  [21], 

for  A0  =  0,  Ho  =  1  and  Go  =  x  (39) 

for  At  0,  Hk  =s  cos(Afci)/ sin( Afc.4)  and  Gk  =  sin(Akx)/sin(AkA)  (40) 


v  — 


(l-/)2 


1  +  / 


The  Ait's  are  complex  parameters  arising  from  this  classical  separation-of-variables  tech¬ 
nique  and  will  be  evaluated  so  as  to  help  satisfy  the  boundary  conditions.  (The  values  so 
determined  will  be  referred  to  as  eigenvalues.)  The  coefficients  of  Ho  and  Go  in  Eqs.( 33-36) 
are  special  perturbation  functions  for  Ao  =  0.  In  the  low-viscosity  theories  of  References 
[4,  10-12]: 


U’o  =  ii0  =  t30  =  po  =  0 

Eqs.(33-36)  can  be  substituted  in  Eqs. (29-32)  to  obtain 

ruk  +  ufc  -  ivk  +  Afclrtek 
(/?e)_l[uk  +  u'Jr  ~(B  +  r'2),}*]  +  2(1  +  x(r2Rc)~l ]0k 
(Re)-'  [uk  +  v i/r  -  (D  +  r“J)ut]  -  2(1  +  t(raJle)~l]tU 

(JZe)-,[u'k  +  w[/r  -  Bwk] 


(41) 

0 

(42) 

p* 

(43) 

-tpk/r 

(44) 

—  XklPk 

(45) 

where 


B  =  r-a  +  Aj-«(l  -  f)Re 


Aol  =•  0 

A  03  =  1 

Akl  =  Xki  —  \k,  k  0 

(  y  =  d(  )/ dr 


These  equations  can  be  put  in  canonical  form  by  the  introduction  of  six  new  sets  of 
variables  Znk  that  satisfy  the  differential  equations  and  a  set  of  constant  multipliers  c*  that 


will  be  used  to  satisfy  die  boundary  conditions 

fekZ\k(r )  = 
fckZik(r)  = 
fckZsk(r)  = 
fckZAk(r )  = 

fckZik(r)  = 

fckZnk(r)  = 


(46) 

tile  -  iVk 

(47) 

K 

(48) 

Wk 

(49) 

wk 

(50) 

Pk 

(51) 

For  these  new  variables,  Eqs. (42-45)  become 


z;*  = 

—  Zlklr  —  AjfciZu 

(52) 

%k  = 

—Z2  k/r  —  A  HZ**  —  iZzk 

(53) 

"3  k  ~ 

2 (Rz  +  ir~2)Zne  +  i(B  4-  r  2)(Z2k  ~  Z\k) 

-Zzkif  —  ( iRe/r)Z6k 

(54) 

Z<k  = 

Zsk 

(55) 

zU  = 

BZAk  —  Z$k/r  —  ^kiR^Zek 

(56) 

z'6k  = 

(Re)-'[-BZlk  +  (2 Me  -  r~2)(Z2k  -  Zlk) 
+iZzk/r  —  AkiZsij 

(57) 

III.  Boundary  Conditions 

The  boundary  condition  on  each  wall  of  the  container  is  that  the  liquid  must  have  the 
same  velocity  as  the  wall: 


On  the  cylindrical  wall  (f  =  1) 


Ui(l,x)  =  *(1  —  /) 

(58) 

u(l,x)  =  -t(l  -  f)x 

(59) 

=  -(1  -  f)x 

(60) 

while  on  the  two  end- walls  (i  =  i  A) 

m(r,±>i)  *  »(l-/)r 

(61) 

(62) 

Si(r,  ±^4)  =  ^(1-/)A 

(63) 

For  a  fully-filled  container,  the  conditions  on  the  axis  are^10' 

lii(0,  x)  m  u(0,  x)  -  tu(0,  x)  =  g(0,  x)  = 

=  0 

(64) 

On  the  free  boundary  of  a  partially  filled  container,  the  tangential  shear  should  vanish  and 

the  normal  stress  should  be  the  inner  pressure,  po. 

(Jk)"1  f^0’1)  %(r«*x)  ~  *u(ro,x)]' 

dr  ro 

=  0 

(65) 

[^(ro,*)  ,  5u(r0,x)' 
{Re)  dr  +  ftr  j 

=  0 

(66) 

£(r„,x)  +  ,r»f(r»'l)-2(Re)-'a“(2'i: 

1  4“  /  UT 

-  =  0 

(67) 

The  no-slip  conditions  at  the  rod  reduce  to: 

Ui(r0,x)  =s  i(l-/)r0 

(68) 

ll(r0,x)  =  -»(l-/)x 

(69) 

li(r0,x)  =  -(l-/)x 

(70) 
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IV.  Low- Viscosity  Solutions 


For  high  Reynolds  numbers  (Re  >  5000),  three  distinct  solutions  have  been  obtained 
to  the  preceding  differential  equations  and  boundary  conditions: 

1.  Stewartson  (S)  inviscid  solution; 

2.  Stewartson- Wedemeyer  (SW)  inviscid  solution  modified  by  wall  boundary  layers; 

3.  Kitchens- Gerber-Sedney  (KGS)  Unear  Navier-Stokes  solution  with  end- wadi  boundary 
layers. 

In  all  three  methods,  Co  =  0  and  the  A0  perturbation  function  for  m  is  i(l  —  /)r,  a  choice 
that  satisfies  the  inhomogeneous  normal  end-wall  boundary  condition,  Eq.(61). 


1.  Stewartson  (S)  Solution 

For  infinite  Reynolds  number,  Eqs. (52-57)  reduce  to: 

Z[k  =  —Z?k/T~^kZ*k 

Z'isk  =  *(1  +  f)Z\k  +  2i Zik 


(1  +  f)Zlu 
1  -/ 
i^kZ«k 

f»(l  +  f)Zlk 


(1  ~f)r 


At 

(l-/)3 


(3  -  /)(1  +  f)Z\k  -  ^ 


The  normal  velocity  at  the  flat  end-walls,  Eq.(61),  cam  be  satisfied  by 

c*  =  0,  Jb  =  2,4,6,...  (77) 

A,  =  JJ.  *-1 .2.3. -N  (78) 

Since  the  perturbation  functions  vanish  for  even  k,  we  wiU  take  N  to  be  odd  and  consider 
only  the  (N+l)/2  odd  k  values  1,3,5,. ..,N  both  in  this  section  and  in  the  discussions  of 
the  SW  and  KGS  solutions  that  follow.  It  is  interesting  to  note  that  the  second-order 
differential  equation  for  that  can  be  obtained  from  Eqs. (71-74): 

r’Z6",  +  tZ'm  +  |  'Ijj'-  -  r>  -  1  j  Z„  =  0  (79) 

is  essentially  Bessel’s  equation  of  order  1.  Thus  the  six  Z:*’s  are  linear  combinations  of  the 
Bessel  functions  J\  and  Y\  and  their  derivatives;  these  combinations  are  given  in  Table  1. 


8 


Next,  the  normal  velocity  at  the  cylindrical  wall,  Eq.(59),  can  be  satisfied  by  Eq.(34) 


when 


where 


=  1,  ck  =  fak 

r  -2*(1  -  f) 


1  +  / 

and  where  the  ajt’s  are  the  coefficients  of  a  least-squares  fit  of  x  to  a  sine  series: 

N 

X  =  Y,akGk(£) 
k=l 

Since  a*  =  0  for  k  even,  the  Cfc’s  of  Eqs.(77)  and  (80)  are  compatible. 

The  second  boundary  condition  for  Eqs. (71-74)  is 
(a)  for  a  fully-filled  cylinder: 


Zek(0)  =  0 


(b)  for  a  partially-filled  cylinder: 


„  ,  x  ,  ir0Zlk(r0)  ir0f 
Z«(r„)+  w 

(c)  for  a  cylinder  with  a  central  rod: 

Zxk(rQ)  =  1 


2.  Stewartson-Wedemeyer  (SW)  Solution 

The  Stewartson  inviscid  solution  couldn’t  satisfy  the  tangential  velocity  conditions  on 
the  container  walls.  In  1965  Wedemeyer  introduced  three  oscillatory  boundary  layers  on 
these  walls.  These  boundary  layers  created  small  viscous  contributions  to  the  liquid  mo¬ 
ment,  but  more  importantly  they  modified  the  larger  inviscid  pressure  moment  term.  The 
boundary  layers  caused  an  outflow  velocity  that  changed  the  normal  velocity  conditions 
for  the  inviscid  solution.  The  presence  of  the  boundary  layers  on  the  end  walls  modifies 
the  eigenvalues,  A*,  and  the  Fourier  coefficients,  ci,  (see  Appendix  C).  The  boundary  layer 
on  the  lateral,  cylindrical  wall  modifies  the  boundary  condition  on  Zu(  1)  (see  Appendix 
D).  Equations  (59)  and  (61)  for  the  inviscid  velocity  become: 


/■,  \  ^1^(1, 3?)  •/.  4, 


where 


Sa  = 


/2(1  —  f)Re 


C  _  -(1+0  l  +  f  ,  »(3 -/)' 

c  "  2(1  -f)y/2Re  [V^1!  ^  v/rTT. 

1  —  en 

e  _  - L  n  =  1  for  a  central  rod ;  2  otherwise 

1  +  e? 

d  =  exp[(r0  -  l)/<5a] 


For  the  central  rod,  a  fourth  boundary  layer  modifies  Eq.(69): 


ii,(ro,-'-)  + 


/  rQ6ae2  \  dll, (r0,  x) 
\ro  +  /  <9r 


-Oi  -  /)* 


(86) 


Eqs. (84-86)  change  Eqs.  (7S,S0,83)  to  the  following: 

cos(AfcA)  +  A k6c  sin(\kA)  —  0,  (k  odd) 
SaClZ'J  1) 


1  - 


=  1 


Z\k(ro)  + 


roSae-2 


Z[k(ro)  =  1 


To  +  6ae7j 

A  good  approximate  solution  of  Eq.(87)  for  the  eigenvalues  is  given  by 

A* 


kit 


2(A-6e) 


(87) 

(88) 

(89) 


(90) 


Eq.(901  yields  estimates  for  solving  Eq.(S7)  iteratively  for  exact  values  of  At.  The  ap¬ 
proximation  (90),  however,  is  only  valid  for  |6C|  <<  A,  in  which  case  the  cosine  in  (87)  is 
nearly  zero.  For  low  Reynolds  numbers,  a  solution  of  Eq.(87)  exists  for  which  the  cosine 
has  significantly  larger  magnitude.!25'  An  adequate  first  estimate  of  this  special  eigenvalue, 
is  given  by 

A,  si  ±i/6c  for  |A,  -  1|  >  0.1  (91) 

where  the  sign  is  chosen  so  that  the  real  part  of  A,  is  positive  (usually  the  minus  sign  for 
r  <  1  and  the  plus  sign  for  r  >  1).  This  special  eigenvalue  and  its  associated  eigenfunctions 
should  be  included  in  all  summations  of  the  eigenfunctions  for  boundary  conditions  and 
moment  coefficients. 


The  SW  solution  gives  excellent  agreement  with  experiments  for  Re  >  5,000. 


3.  Kitchens-Gerber-Sedney  (KGS)  Solution 

The  KGS  solution  employs  the  full  linearized  Navicr-Stokcs  equations  [Eqs. (52-57)], 
but  retains  the  Wedemeyer  end-wall  boundary  layers.  This  allows  the  use  of  the  easily 
computed  SW  eigenvalues  of  Eq.(87).  The  cylindrical  boundary  layer  is  eliminated  and 
conditions  (58-60)  are  satisfied  exactly.  The  boundary  conditions  on  the  Z]k's  at  the  outer 
wall  are: 

Z<k( 1)  =  0  (92) 
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Zlk(  1)  =  1,  ck  =  fak  (93) 

Z2k(  1)  =  0  (94) 

Note  that  ck/ak  is  the  same  as  in  both  the  S  and  SW  solutions.  The  inner  boundary 
conditions  are: 

(a)  for  a  fully-filled  cylinder:2 

Zik(  0)  =  Z<k(  0)  =  Z6k(  0)  =  0 

(b)  for  a  partially-filled  cylinder: 

ro^3t(ro)  —  iZik(ro) 

A*Zii(ro)  -  Zsk(ro) 

Z«,(r„)  +  -  (J-)  [ij3»(r„)  -  A,Z«(r„)j 

(c)  for  a  cylinder  with  a  central  rod: 

Zu(r0)  =  1 
Zik(ro)  —  Z<*(r  o)  =  0 

The  KGS  solution  is  probably  better  than  the  SW  solution  since  it  depends  on  fewer 
boundary-layer  approximations;  however,  KGS  has  been  limited  to  applications  for  which 
Re  >  1000.  An  important  characteristic  of  the  KGS  method  is  the  fact  that  it  can  be 
modified  to  give  an  exact  low- Reyr olds- number  solution  (the  HSG  solution). 


=  0 
=  At 


V.  High-Viscosity  (HSG)  Solution 


In  the  low  viscosity  solutions,  the  k  =  0  solution  was  selected  to  satisfy  the  boundary 
conditions  at  the  flat  end- walls.  The  A*’s  were  then  the  eigenvalues  of  the  differential 
equations  for  functions  of  x.  The  boundary  conditions  on  the  cylindrical  wall  were  then 
satisfied  by  a  Fourier  series  fit  of  these  functions. 

The  HSG  method  reverses  this  process.  The  k  =  0  solution  is  selected  to  satisfy  the 
boundary  conditions  on  the  cylindrical  surface.  The  A*’s  become  the  eigenvalues  of  the 
differential  equations  for  functions  of  r  and  the  end- wall  boundary  conditions  are  satisfied 
by  a  least-squares  fit  t~  eigenfunctions  of  r. 

The  conditions  for  the  k  =  0  perturbation  functions  are: 

(a)  in  general: 


_  Zio(l)  =  1 

3  For  computational  iwm,  condition*  an  not  specified  at  the  singular  point  r  •  0  but  at  some  small  (vw'liv*  value  r®,  Tho 
derivative*  and  can  be  obtained  from  Eq*.(53,S5,57)  to  improve  the  approximation*  used.  Actually,  Ref*r*sc*  10 

used  a  power  series  in  r  with  a*  many  a*  thirty  term*  to  relate  condition*  at  r  m  0  to  thov*  at  r  *  rg.  If  ro  1*  taken  small 
enough,  on*  or  two  terms  Mem  to  be  sufficient. 
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Z2o(l)  —  -^4o(l)  —  0 

CO  =  / 


(b)  for  a  fully-filled  cylinder: 

Z20(  0)  =  Z4  o(0)  =  Zeo(0)  =  0 

(c)  for  a  partially-filled  cylinder: 

r0Z3o(ro)  -  iZ2o(r0)  =  0 
Zio(ro)  —  Z$  o(ro)  —  1 
rr  /  \  ,  ir0Zio(r0)  2iZjo(r0)  irQf 

Zx(r0)  +  —r-j  — - — 

(d)  for  a  cylinder  with  a  central  rod: 

Zio(r0)  =  1 

Zio(ro)  =  Z4o(ro)  —  0 

The  conditions  for  the  k  ^  0  perturbation  functions  are: 

(a)  in  general: 


Zu(  1)  =  Z»(  1)  =  Z4*(  1)  =  0  (95) 

(b)  for  a  fully-filled  cylinder: 

Z2*(0)  =  Z4*(0)  =  Z**(  0)  =  0  (96) 

(c)  for  a  partially-filled  cylinder: 

T-oZ3i(r0)  -  iZik(ro)  =  0  (97) 

AfcZj*(r0)  —  Z?fc(r0)  =  0  (98) 

Z«*(r0)  +  ,r°Vr°).  -  (j|)  [.Z3*(r0)  -  AfcZ,*(r0)]  =  0  (99) 

(d)  for  a  cylinder  with  a  central  rod: 

Z\k(r  o)  =  Zj*(ro)  =  Z4k(ro)  —  0  (100) 


For  add  three  inner  boundary  conditions,  the  solutions  to  Eqs.(52-57)  must  satisfy  the 
homogeneous  boundary  conditions;  hence  these  solutions  are  usually  identically  zero.  We 
can,  however,  select  A*  so  that  other  solutions  axe  possible.  In  particular,  we  will  replace 
the  conditions  Z\k(  1)  =  0  by  Z**(  1)  =  1  and  obtain  a  set  of  solutions  for  particular  values  of 
A*.  This  cam  be  done  by  a  numerical  orthonorrnalization  process  available  in  a  very  generad 
and  flexible  code  called  SUPORTJ26!  This  code  adlows  the  calculation  of  Zu(l)  as  a  function 
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of  Afc.  The  value  of  A*  that  makes  Zu(l)  zero  is  an  eigenvalue  whose  eigenfunctions  satisfy 
all  the  boundary  conditions  on  the  inner  and  outer  cylindrical  surfaces,  r  =  r0, 1. 

The  actual  determination  of  the  eigenfunctions  is  a  delicate  process.  A  simple  Newton 
iterative  method  (strictly  speaking,  the  rtgula  falsi  method^)  has  been  coded;  it  will 
converge  rapidly  to  an  eigenvalue  if  a  good  initial  guess  is  available.  Hall  et  al.  showed^ 
that  for  a  fully-filled  cylinder,  each  eigenvalue  falls  into  one  of  three  families.  Figure  1 
(from  Reference  [21])  shows  72  of  these  eigenvalues  for  Re  =  1000,  r  =  0.1.  (Note  that 
eigenvalue  (3,1)  is  a  maverick  in  this  three-family  concept.) 

For  Re  =  10  and  r  =  0,  the  72  eigenvalues  nearest  the  origin  in  the  half-plane  i?{Afc}  > 
0  were  computed  by  a  numerical  search  process;  these  72  eigenvalues  are  tabulated  in  Table 
2.  From  these  values,  additional  sets  of  eigenvalues  were  generated  by  taking  small  steps 
in  both  Re  and  r.  The  resulting  eigenvalue  data  base  has  been  stored  on  our  VAX/8600 
computer.  For  the  fully- filled  cylinder,  this  base  consists  of  6216  files,  each  file  containing 
72  eigenvalues  for  a  specified  pair  of  values  (Re,r).  There  are  14S  specified  Re  values: 


1 

to 

20 

in  steps  of 

1 

[  20  Re  values] 

30 

to 

100 

in  steps  of 

10 

[  8  Re  values] 

120 

to 

2500 

in  steps  of 

20 

[  120  Re  values] 

and  42  r  values: 

-1 

to 

-.40 

in  steps  of 

.10 

[  7 

values] 

-.35 

to 

-.10 

in  steps  of 

.05 

C  6 

values] 

0 

to 

.10 

in  steps  of 

.01 

C  11 

values] 

.15 

to 

.45 

in  steps  of 

.05 

C  7 

values] 

.50 

to 

1.50 

in  steps  of 

.10 

£  11 

values] 

yielding  the 

148 

II 

C! 

Tf 

X 

6216  files.  By  interpolating 

entries  in  1 

obtains  good  estimates  of  the  72  eigenvalues  when  1  <  Re  <  2500  and  — 1  <  r  <  1.5. 
These  starting  values  allow  computation  of  the  72  eigenvalues  in  less  than  five  minutes  on 
the  VAX. 


Similar  though  not  as  extensive  data  bases  have  been  established  for  the  partially- 
filled  cylinder  and  for  the  cylinder  with  a  central  rod.  For  both  the  partially-filled  and 
central- rod  cases,  files  have  been  created  for  1  <  ile  <  2500  and  0  <  r  <  0.14  and  for 
r0  =  0.1  and  0.3.  Three-way  interpolation  (for  Re,  r  and  ro)  then  yields  good  eigenvalue 
guesses  within  the  specified  parameter  region. 

Once  the  eigenfunctions  are  available,  they  can  be  used  to  satisfy  the  conditions  on 
the  end-walls  by  a  proper  selection  of  the  c*’s: 

N 

]C ckZ\k{r)gk  =  Ri  =  A/[l  -  Zio(r)]  (101) 

fcml 
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.V 

Y^,ckZ2k(r)gk 

k=  1 


.V 

Y^,ckZ*k{r)hk 

/k=l 


/?2  =  -.4/Z»(r) 
R3  =  — /Z.’o(r) 


(102) 

(103) 


where 


^  =  Gt(A) 

/ik  =  ff*M) 

As  can  be  seen  from  Table  2,  the  imaginary  part  of  At  can  be  quite  large  and  the  use  of 
definition  (38)  for  Gk  and  Hk  would  cause  gk  and  hi,  to  have  extremely  large  magnitudes. 
Definition  (40)  avoids  this  difficulty. 

The  Cfc’s  are  selected3  to  give  the  least-squares  fit  of  Eqs.(lGl-103),  where  the  function 
to  be  minimized  is 

R=  ['{ IE  ckZlk  -  Rx  I2  +  1  £  -  R2  |J  +  |  Y.  c^hk  -  R*  IV  dr  (104) 

Jr0 

The  quality  of  the  resulting  fit  is  measured  by  the  error  e,  where 

t  =  (R/Ro)1'2  (105) 

and  where  Ro  is  the  value  of  R  for  c*  =  0: 

R*=  /‘{I  Rx  |2  +  |  |J  4-  !  R3  l2}r  dr 

Jr0 

If  e  <  .05,  the  fit  is  reasonably  good  and  we  have  enough  eigenfunctions  (that  is,  the  value 
assigned  to  N  is  adequate). 


VI.  Equations  for  Two  Immiscible  Liquids 


If  two  immiscible  liquids  are  fully  spun-up  in  a  cylindrical  container  at  zero  coning 
angle,  they  will  be  separated  by  a  cylindrical  interface.  At  this  interface,  r  =  jq ,  the  density 
and  kinematic  viscosity  will  jump  from  the  inner-liquid  values  (p2)u2)  to  the  outer-liquid 
values  (/>i ,  i>\ ).  The  density  px  =  pi  will  be  used  to  nondimensionalize  the  pressure.  The 
symbol  Re  will  denote  the  Reynolds  number  of  the  outer  liquid. 

For  small- amplitude  coning  motion,  the  interface  is  described  by: 


r,  =  r,[l  -  R{rix{i)K'}} 


(106) 


where 


2,(x) 


«u(n,g) 

d-/)r. 


■’Other  pouihle  •election  criteria  for  tha  c» '•  are  given  in  Reference  21.  In  moM  caeca  where  the  fit  ie  good,  the  reeulting 
computed  aide  moment  ie  eeeentiaily  the  tame. 
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The  pressure  relation  (28)  assumes  two  forms,  one  for  each  side  of  the  interface: 


P~Po 

pia24>7 


H  =  9.  0  (or  r,  or  r,)  <  r  <  « 

P  L  Z 


=  ft 


r?  —  Tnl  r3  —  r? 


+  _  7J{£A'*},  r,  <  r  <  rc 


where 


(107) 


ft  =  pi(p\ 

Differential  equations  (52-57)  become 

Z\k  ~  ~Z‘iklr  —  ^k\^4k 
Z'2k  =  —Zik!r  ~  ^lel  Z\k  ~  *Z3 k 

Z'3k  =  2 (qirRt  +  ir~2)Z\k  +  *(i?  +  r  2)(Zi k  ~  Zit) 

-Z3*/r  -  (iqirRc/r)Z6k 

Z\k  ~  Z>k 

Z'ik  =  BZ41,  —  Zsk/r  —  XkiqirReZak 

Z'sk  =  (?2rJRerl(~5Zlk  +  (2i9:r/2e-r-2)(7J*-Zu) 

+*Z3k/r  —  XkiZsk] 

where 

5  =  r-3  +  Aj  —  t(l  -  f)qirRe 

qir  =*  «/a  5  ^/i/j  r0  <  r  <  r, 

=  1,  rj  <  r 

At  the  interface,  the  velocities  and  stresses  are  continuous;  therefore, 

IBK,*)  =  Ui(r*,x) 

lt(r,",x)  =  n(r+,x) 

li(r,",x)  =  li(r+,j) 

/«  v-i  r^(ri',x)  »[li(rf,x)-ic(rr,*)]l  f^(rf,x)  t[ii(r+,x)- 

*(ne)  l-* - n - j  =  (i!e)  [— 5 - n 

.  fSutOf,*)  $tt(rf,x)l  .  dui(rf , x)  ,  9tt(rf,x)' 

93(/?e)  3;:  +  ~qx  =  ^  [— —+~"dr~\ 


(108) 

(109) 

(110) 
(111) 
(112) 

(113) 


(114) 

(115) 

(116) 
»'tt(ri+,x)l 

(117) 

(118) 


«  [jKrr.xi  +  fia +!a^_2(/fc)-.MjLi> 

(119) 

where 

„  _  ft  _  Pa 

ft  ~  1 

ft  Pi 

the  ratio  of  th*.  inner  and  outer  dynamic  viscosities.  Thus  the  pressure  perturbation  g  and 
the  radial  derivatives  of  ii  and  ui  have  jump  discontinuities  at  the  interface. 
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VII.  Binary  Solutions 


1.  Stewartson  Inviscid  Binary  (SB)  Solution 


Eqs. (71-83)  still  apply,  with  the  additional  complication  of  a  jump  in  the  pressure  term 
at  the  interface  r  =  rx.  If  we  denote  the  binary  perturbation  functions  for  the  pressure 
and  velocity  bv  Z\kB  and  Ze*s,  Eqs. (115,119)  reduce  to 


where 


&ZekB 


0 

(<?.-!) 


ZekB^x  )  + 


irxZuB(rs  )  irj 

W  2  . 


A(  )  =  (  )T=rt  -  (  )r„r 


(120) 

(121) 


These  binary  perturbation  functions  can  be  constructed  from  the  single-liquid  func¬ 
tions  by  the  following  device: 


Z'us(r)  =  Zik(r)  +  akZrk(r), 

r  <ri 

—  Zu,(i')  +  bkZri(r), 

r  >  r. 

(122) 

Z6kfl(r)  =  Z6k(r)  +  akZui(r), 

r  <  rt 

=  Z6*(r)  +  SkZ,2k(r), 

r  >  rt 

(123) 

where  the  Zrk’s  and  Z12*’s  are  solutions  of  Eqs. (71-74)  when  each  subscript  ‘nk’  has  been 
replaced  by  4(n+6)k’  and  the  boundary  conditions  are 

(a)  in  general: 

^7*(fi)  =  1 

ZTk(l)  =  0 


(b)  for  a  fully-filled  cylinder: 

Zi7k(  0)  =  o 


(c)  for  a  partially-filled  cylinder: 

r,»(r, o)  +  'IsMpl  =  o 

(d)  for  a  cylinder  with  a  central  rod: 

Zt  k(r0)  =  0 

Note  that  Z12k  has  a  jump  discontinuity  at  The  Z7*  and  Zm  functions  can  be 
computed  by  two  calls  to  the  SUPORT  code:  for  0  <  r  <  r,  (or  r0  <  r  <  rt)  and  for 
<  r  <  1.  The  ak’s  and  Sk’s  are  determined  by  Eqs. (120-121): 

a*  =  bk  (124) 
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N 


k 


Z\ )  -  qiZMr; )  + 


»(1  ~  9i)ri 

W 


=  (?«-!) 


Z$k{r\)  + 


iriZikln) 

W 


(125) 


Equations  •imilar  to  (124-125)  were  used  by  Scott!28'  to  obtain  eigenfrequencies  r*  for  which 
the  inviscid  liquid  moment  was  infinite. 


2.  Stewartson-Wedemeyer  Binary  (SWB)  Solution 


If  //,  ^  u\,  the  SWB  method  must  use  an  average  kinematic  viscosity;  this  average 

was  chosen  to  be  „ 

^(r?-r2)  +  ^(l-r2) 


v*  = 


1  -  To 


(126) 


From  this,  an  average  Reynolds  number,  7?c„,  and  an  average  value,  6n,  of  Se  can  be 
determined;  the  eigenvalues  can  then  be  computed  from  Eq.(S7)  for  this  average  value.  An 
approximation  for  these  eigenvalues  is 

<127> 


Next,  two  additional  boundary  layers  are  inserted  at  the  interface  to  satisfy  the  continuity 
of  the  tangential  stress,  Eqs.(65-G6).  This,  however,  requires  a  discontinuity  in  the  inviscid 
radial  velocity,  since  the  outflow  in  the  two  boundary  layers  are  not  equal.  Eqs. (120-121) 
have  new  jump  values  in  terms  of  the  Stewartson  pressure  jump  (A Z^kB)s- 


where 


A  Z\kB  =  (1  —  <Ji)&*F 
AZats  »  (A 


F  = 


e3  * 
= 

9s  = 


{ 


U  +  ^)Az;t- 


2r3 

\d(rZik) 

_  2t£3  1 

ri 

dr 

r»  / 

exp((r,  -  1)/*,) 

+  9»  +  (1  -  9s)e3]  -  (1  -  9t)(l  +  «3)*« 
9i/92/3 


(128) 

(129) 


Note  that  Eqs.(128-129)  reduce  to  Eqs.(120-121)  for  Sa  =  0.  The  SWB  equations  in 
this  report  are  based  on  the  assumption  that  the  radial  thickness  of  the  inner  liquid  is 
substantially  greater  than  the  boundary  layer  thickness;  that  is,  r,  -  r0  >>  |5„|.  Reference 
7  considers  the  case  where  this  assumption  is  invalid. 


Finally,  the  presence  of  boundary  layers  on  the  cylinder  wall  and  on  the  rod  requires 
new  conditions  on  Z\k(  1)  and  Z7*(l)  for  all  cases: 


Z,‘(1)“(r^l)Z;‘(1)  =  1  -  (130) 

£;,(!)  —  ^  j  Z'7k(l)  =  -(2c.v/jia/f,)A (131) 
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where 


c.s  =  [1  +  ?s  -  (1  -  7')f3]/[l  +  9s  +  (1  ~  flska) 
Cti  =  [1  +  75  +  (i  —  qi)cl\[l  —  f5^a] 
and  on  Zu-(ro)  and  Zri(  r0)  for  the  rod: 

z"{r°)+ (^rk)  = 1 

2,*<ro)+(^:)z;‘(r“)=0 


(132) 

(133) 


where  St  —  SaqjI/:l. 

Murphy^  used  equations  similar  to  (128)  and  (129)  to  compute  liquid  moments  that 
showed  excellent  agreement  with  experiment. 


3.  KGS  Binary  (KGSB)  Solution 

The  KGSB  solution  uses  the  simple  eigenvalues  of  Eq.(127)  but  requires  the  use  of 
jump  conditions  (117-118)  in  addition  to  the  pressure  jump  conditions  given  by  Eq.(119) 
or  (121).  The  three  jump  conditions  for  k  =  1,2 . N  are 

Z3tB(r1+)  -  Z3j,a(r,~)  =  A ZzkB 
Za,B(rt )  —  Zs*s(r  )  =  A  Zs*a 
Z6*a( rf)-  Z6*fl(rC)  =  AZ^a 

where 

A Z3*a  =  (73  -  l)(Z3*a(cC )  —  (ZjtaCfi )/ri] 

A Z5*a  =  (73  -  1  )[Z5*:s( r*i" )  “  )  +  At] 

and  where  AZg*a  is  defined  in  Eq. ( 121). 

In  order  to  satisfy  conditions  (134-136),  we  need  to  introduce  three  sets  of  six  auxiliary 
functions,  defined  in  Table  3.  These  functions  are  continuous  in  the  velocity  components. 
Substitution  in  Eqs.(134-136)  yields  the  values  of  <4,  e*,  and  fk.  The  resulting  values  of 
(Z,jta)i=i...fl  can  then  be  used  to  compute  the  liquid  moment  exerted  by  a  binary  liquid 
payload. 


(  34) 
( *35) 
■  136) 


4.  HSG  Binary  (HSGB)  Solution 

In  the  KGSB  solution,  wo,uo,0o,  and  po  were  zero  and  the  auxiliary  functions  of 
Table  3  were  used  for  positive  k.  For  HSGB,  the  binary  A<>  eigenfunctions  must  also  be 
computed  by  use  of  the  auxiliary  functions  given  in  Table  3.  The  jump  conditions  for  the 
A0  eigenfunctions  are 

AZioa  =  (73  -  l)(-^3oa('”j' )  -  iZ2oa(ri )/rt]  (137) 


18 


AZsos  =  (<73  —  l)[Z5os(rf )  —  Zioa(ri)  +  1] 

AZeoa  —  (<7i  —  1) 


„  ,  _N  .  ir\ZWB{rx)  ,  irxf 

^6oa(r,  )  +  — — j —  +  — 


2(g3  —  l)^{oa(rO 
Re 


The  jump  conditions  for  the  other  eigenfunctions  (k  =  1,2,...,  N)  are 

A Zmb  =  (<73  -  l)(^3*a(rj’)  —  *^2*a(ri )/ri ] 

AZskB  =  (?3  -  l)[-^5fca(»*r )  ~  ^fcZi*a(ri)] 


A  Z6*a  =  (<7i  -  1) 


Zataifj')  f 


»>i^ua(ri) 


1-/ 


2(g3  ~  l)Z'(fca(r1) 
Re 


(138) 

(139) 


(140) 

(141) 

(142) 


These  jump  conditions  are  used  to  determine  the  necessary  dt’s,  e*’s  and  /*’ s.  The 
resulting  binary  eigenfunctions  are  then  used  to  satisfy  the  end-wall  boundary  conditions 
of  Eqs.(101-103);  that  is, 


JV 


53  ckZ\kB(r)gk 

km  1 

=  Rx  =  Af[l  -  Z10a(r)] 

(143) 

N 

53  ckZjkB(r)gk 

leal 

=  R7  =  - AfZ20B(r ) 

(144) 

53  ckZ*kB{r)hk 

km  l 

=  i?3  =  -/Zeoa(r) 

(145) 

The  Ck'a  are  then  selected  to  fit  Eqs.(143-145)  so  that  R  is  minimized,  where 

R~jT9  53c*^l*0  -  R\  |J  +  |  53c*^J*s  -  Ri  |2  +  |  5Zc*Z<ia/»t  —  R3  j2}r  dr  (146) 
and  where 


<l\r  =  qi  (r  <  r,) 

=  1  (r  >  r,) 

The  error  of  the  fit  is  computed  by  Eq.(105)  with 

=  /  *i,{|  Ri  |3  +  j  R,  |J  +  |  R3  |3}r  dr 


VIII.  Liquid  Moment 


If  a  projectile  is  performing  simple  coning  motion,  the  moment  exerted  on  it  by  a  fully 
spun-up  liquid  payload  can  be  defined  in  terms  of  three  dimensionless  moment  coefficients. 
These  coefficients  are  related  to  the  liquid  moment  (A/*e,  A/ye,  Mte)  in  the  coning  coordinate 
system  as  follows: 


A/re  =  A’3 

A/ye  +  iA/,e  =  rnL(^1<^1‘r\CisSt  +  iCust\I\ 


(147) 

(148) 


where 


mi  —  2n  a7cpL 

These  three  coefficients  -  the  liquid  roll  moment  coefficient  Clrm,  the  liquid  side  moment 
coefficient  Clsm,  and  the  liquid  in-plane  moment  coefficient  Cum  -  are  all  functions  of  A, 
r,  e,  Re  and  A'2.  The  liquid  side  moment  affects  the  damping  of  the  coning  motion  while 
the  liquid  in-plane  moment  can  cause  a  frequency  shift.  Thus  our  primary  interest  is  in  the 
liquid  side  moment  coefficient  Clsm  since  it  determines  any  instabilities  that  could  occur. 


The  linear  pressure  at  any  point  of  the  cylinder  can  be  expressed  in  the  coning  coor¬ 
dinates  by  combining  Eqs.(10,ll,107): 


P~Po 

pLa?4>2 


R{(Z+xr)IC}  , 


0  (or  f/  or  fr)  <  r  <  f, 


R{( 2  +  xf)A'*},  n  <  r  <  1 


(149) 


The  linear  coefficient  of  A'*  is  defined  as  the  pressure  coefficient  Cp  and  can  be  expressed 
in  terms  of  Z6k{r)  by  use  of  Eqs.(36,51): 


=  *„/[( 2  -  f)xr  +  £  ckZ*k{r)Gk(x)]  (150) 

*=o 

Note  that  the  distinction  between  r,r  and  x,f  vanishes  in  the  linear  terms. 

The  linear  coefficients  can  be  computed  from  the  fluid  mechanics  theory  of  this  report 
by  integrating  the  pressures  and  viscous  stresses  over  the  lateral  cylindrical  wall  and  the 
two  flat  end-walls.  If  a  central  rod  is  present,  its  pressure  and  viscous  contributions  must 
also  be  included.  Using  subscripts  p  for  pressure,  v  for  viscous,  l  for  lateral  wall,  e  for  end 
wall  and  r  for  rod,  we  have^,12l 

Clsm  +  iCuM  —  Cmpt  +  Cmp*  "h  Cmpr  +  Cmv,  +  Cm„  +  (151) 

Expressions  for  the  components  in  Eq.(151)  are  given  in  Table  4  for  HSG  and  HSGB.  The 
table  also  applies  to  the  S  and  SB  theories  when  the  Reynolds  number  is  infinite,  but 
additional  terms  must  be  introduced  for  the  boundary  layer  theories  (SW,  SWB,  KGS, 
KGSB),  as  indicated  in  Table  5. 

The  linear  fluid  moment!  can  also  be  computed  directly  by  differentiating  the  volume 
integral  of  the  angular  momentum  of  all  the  liquid  in  the  cylinder.  To  do  this,  we  first 
express  the  position  and  velocity  vectors  in  earth-fixed  cylindrical  coordinates  (er,ef,c»): 


R 

= 

(xe,  +  rcT)a 

(152) 

V 

31 

Vte,  +  Vrer  +  Vt£, 

(153) 

e\ 

* 

eycos#  +  e,  sin# 

a 

( - 1\  cte  +  7eyc )  cos  ip  +  e,c  sin  tf* 

(154) 

U 

a 

e,  x  ey 

- 

(A'c«  -  -yc^ain^  +  eiecos0 

(155) 
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Then 


(156) 


R  x  V  =  [rVg&x  -  xVgir  +  (xVT  -  rV,)i,]a 
The  angular  momentum  is 

L  =  a3  j  j  J  pR  xVrdfdx  d\j> 

—  L\esc  +  LiCyc  +  LztJC  (157) 

where  expressions  for  the  Lj  axe  given  in  Table  6. 

The  moment  exerted  by  the  liquid  on  the  projectile  is  the  negative  of  the  derivative 
of  the  liquid’s  angular  momentum: 

M  =  -Z 

—  +  i tCye  +  +  Q  x  L)  (158) 

For  steady-state  constant-amplitude  liquid  motion, 

Lj  =  e  =  0  (159) 

M  —  —  r^e*  x  L 

—  r4>[-KLzixe  +  7  Lztyc  +  (KL\  —  7  L-i)ite\  (150) 

Eq.(160)  -  applied  to  definitions  (147-148)  -  yields  the  very  important  result  that 

7  Clrm  =*  —  C'/Os.w  (161) 

This  result  has  been  obtained  elsewhere^9’1*!,  but  the  above  is  the  simplest,  most  direct 
derivation  and  is  based  solely  on  the  assumption  of  steady-state,  constant-amplitude  mo¬ 
tion  given  by  Eq.(159). 

The  first  terms  in  an  expansion  in  K  of  the  Lj'a,  obtained  by  use  of  Eqs.(16-18),  are 
given  in  Table  6.  Substituting  these  small-K  expressions  for  the  Lj' s  in  Eq.(158),  we  obtain 

Clrm  w  -i2{S}  (162) 

Cisu  +  iCuM  =  (1  -  *e)(S  +  2iB\)  (163) 

Since  Clrm  is  a  quadratic  coefficient,  it  can’t  generally  be  computed  from  a  linear  theory. 
Eq.(162)  neglects  the  nonlinear  part  of  L\.  For  L\  =  0,  Eq.(162)  gives  the  correct  quadratic 
coefficient. 

It  is  easy  to  show  that  S  vanishes  for  f  *  0.  Therefore,  for  r  *»  e  m  0, 

Clsm  =*  Cuu4  =*  0  (164) 

Cum  *  2BX  (165) 

Values  of  Clsm  +  iCuM  obtained  by  this  angular  momentum  approach  will  be  denoted  as 
Cmm.  The  appropriate  formulas  are  given  in  Table  7;  the  derivation  of  these  formulas  is 
given  in  Appendix  B. 
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IX.  Numerical  Results 


1.  Surface  Integral  versus  Volume  Integral 

In  the  HSG,  KGS  and  SW  theories,  the  liquid  moments  were  originally  computed 
by  integrating  the  pressure  and  wall  shears  over  the  walls  of  the  container  (see  Tables  4 
and  5).  The  volume  integral  method  (Table  7)  should  give  the  same  results  if  the  velocity 
perturbation  functions  satisfy  the  linearized  Navier-Stokes  equations  and  the  boundary 
conditions.  For  a  sufficient  number  of  eigenfunctions,  this  is  the  case  for  the  HSG  theory: 
the  two  values  of  Clsm  are  equal.  In  Figure  2,  Clsm  values  obtained  from  HSG  by  both 
the  surface  and  volume  integral  methods  are  plotted  versus  the  number  of  eigenvalues  used, 
for  Re  =  500,  A  =  3.1  and  r  =  0.1.  For  28  eigenvalues  or  more,  the  values  from  the  two 
methods  are  essentially  equal.  For  a  smaller  number  of  eigenvalues,  the  two  values  differ, 
but  the  volume  integral  result  is  always  closer  to  their  common  value  for  more  than  28 
eigenvalues.  This  figure  also  shows  the  error  of  the  fit  to  the  end-wall  boundary  condition; 
we  see  that  the  Clsm  are  the  same  when  the  error  is  less  than  0.1G. 

Gansl29!  uses  the  volume  integral  method  in  his  version  of  SW  and  he  notes  that  this 
method  correctly  predicts  that  the  side  moment  coefficient  goes  to  zero  as  r  goes  to  zero, 
whereas  SW  results  based  on  the  surface  integral  method  do  not  have  this  property.  Both 
Rosenblatt in  his  finite  element  calculations  and  Herbert^9!  in  his  spectral  method  use 
the  volume  integral  approach.  Herbert  emphasizes  the  value  of  the  volume  integral  in 
giving  better  results  for  the  approximate  theories.  Indeed,  for  the  approximate  boundary 
layer  theories  at  low  Reynolds  numbers,  the  average  error  in  velocity  predictions  should 
be  much  less  than  the  error  in  the  velocity  gradients  at  the  boundary. 

In  Figure  3,  the  volume  and  surface  integral  results  for  KGS  and  SW  are  compared  at 
Re  ss  500,  A  =  3  and  r  ranging  from  0  to  0.20.  The  exact  HSG  values  are  also  given  in  the 
figure.  The  surface  integral  results  are  essentially  worthless,  but  the  volume  integral  results 
are  quite  good,  with  KGS  somewhat  better  than  SW.  For  this  reason,  we  will  routinely 
use  the  volume  integral  to  compute  the  liquid  moments  for  all  theories. 


2.  Fully  Filled  Cylinders 

Historically,  our  interest  in  the  moment  produced  by  highly  viscous  liquids  is  based  on 
experiments  by  Miles  Miller!31!.  He  measured  the  liquid  roll  moment  in  a  coning,  spinning, 
liquid-filled  cylinder  and  showed  that  at  a  Reynolds  number  of  about  150,  the  liquid  roll 
moment  was  large  enough  to  cause  flight  instabilities!31!.  The  theoretical  basis  for  this 
very  important  result  is  Eq.(161),  which  directly  relates  the  liquid  roll  moment  of  Miller’s 
experiments  to  the  liquid  side  moment  that  controls  flight  instabilities. 

In  Figure  4,  we  plot  the  negative  of  Miller’s!33-33!^**  versus  logl0Re  for  a  fully 
filled  cylinder  of  fineness  ratio  4.5  and  r  a  0.1.  This  figure  also  gives  the  HSG,  KGS  and 
SW  curves  for  Clsm <  For  HSG  and  KGS,  we  see  good  qualitative  agreement  but  poor 
quantitative  agreement  with  the  Miller  data  for  Re  less  than  about  20. 

The  most  striking  feature  revealed  by  Figure  4  is  the  exceptional  ability  of  the  KGS 
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theory  to  predict  the  much  more  soundly  based  HSG  curve.  The  three  theories  are  com¬ 
pared  for  two  other  fineness  ratios,  3  and  1.5,  in  Figures  5  and  6,  respectively,  and  we 
see  that  the  excellent  agreement  for  the  KGS  theory  persists  for  all  Reynolds  numbers 
greater  than  unity.  The  SW  theory,  however,  does  surprisingly  well  down  to  an  Re  of 
100.  Previous  comparisons  showed  poor  agreement  for  Re  less  than  1000.  The  primary 
difference  between  the  earlier  computations  and  those  of  Figures  4-6  is  that  Figures  4-6 
were  obtained  by  (a)  using  the  volume  integral  method,  (b)  adding  the  special  eigenvalue 
(see  page  10),  and  (c)  imposing  a  more  exact  boundary  condition,  Eq.(Cl9). 

The  SW  method  results  in  a  very  fast  computer  code  since  it  is  based  on  linear 
combinations  of  easily  generated  Bessel  functions.  Indeed,  three-dimensional  surfaces  of 
Clsm  as  functions  of  any  two  of  the  parameters  (r,  e,  Re,  A,  r0)  can  be  obtained  quickly  in 
the  interactive  mode  on  our  VAX  8600. 

The  KGS  method  solves  a  sixth-order  complex  differential  equation  numerically  for 
a  set  of  eigenvalues.  These  eigenvalues  are  quite  easy  to  compute;  usually,  less  than 
twelve  are  required  to  satisfy  the  boundary  conditions  adequately.  For  Re  less  than  1000, 
a  single  calculation  of  Clsm  takes  less  than  15  seconds  on  our  VAX,  but  the  run  time 
increases  with  Reynolds  number.  The  very  exact  HSG  method  unfortunately  requires 
many  more  eigenvalues  (our  code  has  an  arbitrary,  but  easily  extendable,  limit  of  72). 
These  eigenvalues  are  computed  from  tables  of  starting  estimates  and  an  iterative  Newton 
scheme. 

The  agreement  between  the  theories  as  demonstrated  by  Figures  4-6  has  a  very  im¬ 
portant  pay-off.  For  Reynolds  numbers  greater  than  300,  the  improved  SW  method  of  this 
report  can  be  used  to  compute  Clsm  in  a  very  fast  interactive  code.  For  Re  between  0 
and  300,  HSG  results  can  be  reasonably  well  approximated  by  the  KGS  code,  a  process 
somewhat  slower  than  the  SW  program  but  still  fully  interactive. 

This  ability  to  obtain  Clsm  rapidly  over  a  wide  spectrum  of  Re  values  is  demonstrated 
by  Figure  7,  a  three-dimensional  plot  of  Clsm •  versus  Re  and  r,  for  a  fully  filled  cylinder 
of  fineness  ratio  4.5.  The  plot  was  obtained  from  the  KGS  and  SW  codes.  For  the  higher 
Reynolds  numbers,  we  see  the  characteristic  resonance  peak  at  about  r  =  0.20.  For  a  fixed 
t,  a  mild  peak  appears  near  Re  =  50.  Similar  3D  figures  were  first  made  by  Miller*32*. 

An  important  feature  of  Figure  3  that  has  not  yet  been  addressed  is  the  linearity  of 
Clsm  with  r  for  r  less  than  0.2.  According  to  Eq.(148),  a  linear  Clsm  implies  a  side 
moment  that  varies  as  the  square  of  the  coning  frequency  and  is  independent  of  the  spin 
rate. 


The  region  of  linearity  in  r  can  be  characterized  by  the  slopes  ks  and  where 

Clsm  =  ksr,  Clrm  -  k^r 

and  by  rm,  the  maximum  value  of  r  for  which  linearity  is  valid.  We  will  define  rm  to  be 
the  lowest  positive  r  value  for  which  Clsm/t  falls  outside  the  interval  0.9*5  to  1.1*5.  In 
Figures  8  and  9,  we  plot  *5  and  rm,  respectively,  versus  Re  for  Re  between  0  and  1000. 
As  can  be  seen  from  Figure  9,  for  Re  <  250  the  side  moment  coefficient  is  linear  with  t 
for  a  substantial  r  range  (viz.,  r  <  0.2).  This  is  precisely  the  experimental  observation  by 
Miller*31!  for  the  liquid  roll  moment. 
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3.  Partially  Filled  Cylinders  and  Central  Rods 

Stewaxtson’s  original  inviscid  world1!  and  Wedemeyer’s  viscous  extension!2!  covered 
both  fully  and  partially  filled  cylinders.  Frasier!3!  extended  the  SW  theory  to  a  cylindrical 
cavity  with  a  central  rod.  Both  the  KGS  and  HSG  theories!10-13'21!  were  originally  restricted 
to  fully  filled  cylindrical  cavities.  With  the  help  of  the  very  general  SUPORT  subroutine!26!, 
it  is  quite  easy  to  extend  the  KGS  theory  to  partially  filled  cavities  and  central  rods;  this 
has  been  done  for  this  report.  As  mentioned  earlier,  the  extension  of  the  HSG  theory  to 
these  cases  requires  the  forming  of  additional  files  of  eigenvalues.  These  eigenvalues  have 
been  obtained,  using  BRL’s  CRAY  supercomputers. 

For  Figure  10,  the  side  moment  coefficient  has  been  computed  by  ail  three  theories  for 
a  partially  filled  cylinder  of  fineness  ratio  3,  with  a  Reynolds  number  of  500  and  r  =  0.10, 
and  plotted  against  the  inner  cylindrical  free  surface  ratio  b/a.  As  we  would  expect,  the 
side  moment  vanishes  when  this  ratio  is  1.0  and  the  cavity  is  empty.  The  generation  of 
HSG  eigenfunctions  gets  to  be  too  tedious  for  b/a  >  0.9,  but  we  see  that  the  KGS  theory 
is  a  good  approximation  to  the  more  exact  HSG  theory  in  the  range  0  <  b/a  <  0.9. 

Figure  11  is  a  similar  plot  for  the  central  rod  ratio  d/a  varying  from  0  to  1.  The  HSG 
eigenfunctions  were  only  calculated  for  d/a  up  to  0.5.  Once  again,  the  computationally 
simpler  KGS  theory  is  a  good  approximation  to  the  HSG  theory.  Thus  for  most  engineering 
calculations,  the  KGS  theory  would  be  the  one  to  use. 


4.  Binary  Liquids 

The  SW  theory  was  extended  to  binary  liquids  (SWB)  in  Reference  7.  The  extension 
of  the  KGS  and  HSG  theories  to  their  binary  versions,  KGSB  and  HSGB,  has  been  done 
as  part  of  this  report.  The  eigenfunctions  and  boundary  conditions  for  the  KGSB  theory 
are  essentially  the  same  as  those  for  the  KGS  theory;  hence  the  eigenfunctions  have  no 
difficulty  fitting  the  lateral  wall  condition,  Eq.(59).  The  HSG  end-wall  boundary  condition, 
Eq.(104),  can  be  reasonably  well  satisfied  by  ten  HSG  eigenfunctions  at  Re  =  10  and  by 
less  than  72  HSG  eigenfunctions  for  Re  less  than  2000.  Unfortunately,  this  is  not  the  case 
for  the  HSGB  boundary  condition  of  Eq.(146)  and  the  HSGB  eigenfunctions.  Even  for 
quite  low  Reynolds  numbers,  a  large  number  of  eigenfunctions  is  required.  For  this  reason, 
we  have  stopped  (at  least  temporarily)  running  our  HSGB  program. 

The  great  success  of  the  KGS  theory,  however,  allows  us  to  place  primary  reliance 
on  the  KGSB  theory  to  calculate  the  effect  of  binary  liquids  of  different  densities  and/or 
viscosities  on  the  liquid  side  moment.  For  high  Reynolds  number,  good  experimental  agree¬ 
ment  has  been  obtained!3"4!;  the  validity  of  the  KGSB  theory  at  lower  Reynolds  numbers 
will  depend  on  future  experiments. 

In  Figure  12,  Cism  is  plotted  against  r  for  a  cylindrical  cavity  containing  approxi¬ 
mately  equal  amounts  of  liquid  (r2  =  0.49),  the  densities  of  the  two  liquids  differing  by  60% 
( P 2  =  0-4pi)  while  the  kinematic  viscosities  are  equal.  The  two  curves  are  for  two  kinematic 
viscosity  values,  corresponding  to  Reynolds  numbers  of  300  and  500.  The  interface  is  at 
rj  =s  0.7,  the  fineness  ratio  is  3. 
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In  Figure  13,  we  restrict  ourselves  to  Re  =  500,  A  =  3,  and  r  =  0.15,  but  vary  the 
density  ratio  q\  from  0  to  1.  The  two  curves  are  for  rq  =  0.5  and  ri  =  0.7.  The  HSG  value 
for  one  liquid  (^i  =  1)  is  plotted  and  we  see  that  the  agreement  is  quite  good.  (It  is  a  wise 
practice  to  confirm  KGSB  calculations,  where  possible,  with  an  HSG  check  point.) 

Finally,  in  Figure  14  we  consider  a  case  where  the  inner  liquid  is  very  viscous  and 
comprises  98%  of  the  cylinder  volume  (rq  —  0.99).  The  small  amount  of  outer  liquid  has 
the  same  density  but  quite  different  viscosity.  (In  any  experiment,  the  density  of  the  outer 
liquid  would  have  to  be  slightly  greater  than  that  of  the  principal  liquid  so  that  the  outer 
liquid  would  be  located  near  the  cylindrical  wall.)  According  to  Figure  4,  for  A  =  4.5 
and  r  =  0.1,  the  maximum  side  moment  coefficient  occurs  at  an  Re  of  approximately  50. 
Hence  in  Figure  14,  we  set  A  =  4.5,  r  =  0.1  and  Re(inner  liquid)  =  50  and  we  vary  qi  from 
0.1  to  190.  Note  that  the  side  moment  coefficient  essentially  vanishes  when  the  Reynolds 
number  of  the  outer  liquid  is  5000.  This  qualitative  behavior  has  been  observed  by  Miller. 
It  is  interesting  to  note  that  for  the  second  curve  (A  =  3),  the  side  moment  is  relatively 
insensitive  to  q?. 


X.  Summary 

1.  The  four  liquid  payload  theories  -  S,  SW,  KGS  and  HSG  -  have  been  developed 
here  in  a  single  theoretical  framework. 

2.  A  very  simple  and  elegant  derivation  of  the  relation  between  liquid  side  moment 
and  liquid  roll  moment  is  given. 

3.  Moment  coefficients  have  been  computed  by  both  surface  integral  and  volume 
integral  and  the  clear  superiority  of  the  volume  integral  is  demonstrated. 

4.  Eigenvalue  tables  for  a  wide  range  of  r’s  and  Reynolds  numbers  have  been  generated 
to  allow  quicker  HSG  calculation  of  the  side  moment  coefficient.  These  tables  originally 
applied  only  to  a  fully  filled  cylindrical  cavity  but  they  have  been  expanded  to  include 
partially  filled  cylinders  and  cylinders  with  a  central  rod. 

5.  Side  moment  coefficients  computed  from  the  more  convenient  KGS  theory  are 
shown  to  be  good  approximations  to  those  obtained  by  the  more  exact  HSG  theory.  This 
agreement  is  shown  to  exist  not  only  for  fully  filled  cylinders  but  for  partially  filled  cylinders 
and  cylinders  with  a  central  rod  as  well. 

6.  Although  the  binary  version  of  HSG  often  fails  to  satisfy  the  end-wall  boundary 
condition,  the  binary  version  of  KGS  does  not  have  this  difficulty  and  gives  good  agreement 
with  experiments  at  high  Reynolds  numbers. 

7.  KGSB,  the  binary  version  of  KGS,  was  used  to  determine  the  effect  of  a  small 
amount  of  lower  viscosity  liquid  on  the  side  moment.  The  theory  predicts  the  large  decrease 
in  side  moment  that  has  been  observed  experimentally. 
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Table  1.  The  Scewartson  Zlk  s  as  Bessel  Functions 


zlk  =  i^rb[EkJo  +  FkYo]  +  ^ [Ekj\  +  FkYx } 
Zu  =  -'-£f[EkJ0  +  FkY0 }  +  (3 IfclEkJt  +  FkYt] 


Z3k  =  -^r[EkJ0  +  FkYo}  + 


7  _  nl 

(3-/)r*  (l-/)2 


[EtJi  +  F*r,] 


z4*  =  +  iW] 

Z5fc  =  ^[EfcJo  +  nr0]  -  jfzj)'.[EkJi  Hr  FkY{\ 
z6k  =  EkJy  +  FkYx 


where 

/i  =  ^(3 -/)(!  + /)/(!-/) 

•/n  =  Jn(fb^kr) 


=  Bessel  function  of  the  first  kind  of  order  n 

K  =  Yn(fb\kr) 

—  Bessel  function  of  the  second  kind  of  order  n 

and  where  the  F^’s  and  Fit’s  are  determined  by  the  boundary  conditions. 

The  condition  Z\k{  1)  =  1,  for  example,  requires  that 

((1  -  f)fb\kJo«  +  (1  +  f)Ju]Ek  +  [(1  -  f)fb\kY0a  +  (1  +  f)Yla)Fk  =  — 1(3  -  /)( 1  +  /) 

where  the  subscript  ‘a’  on  a  Bessel  function  denotes  the  value  at  r  =  1. 

For  the  fully-filled  case,  the  condition  Zn ic(0)  =  0  requires  that  Fk  =  0. 


Table  2.  Eigenvalues  for  a  fully-filled  cylinder,  Re  =  10,  r 


n 

-  Family  1  - 

-  Family  2  - 

-  Family  3  - 

i 

.904,  4.290 

1.962, 

-3.697 

1.309,  1.518 

2 

.573,  7.528 

1.925, 

-6.702 

2.837,  6.690 

3 

.415,  10.796 

1.983, 

-9.791 

2.567,  9.787 

4 

.326,  14.052 

2.054, 

-12.891 

2.484,  12.895 

5 

.269,  17.289 

2.124, 

-15.997 

2.463,  16.005 

6 

.229,  20.510 

2.1SS, 

-19.109 

2.468,  19.117 

7 

.200,  23.718 

2.247, 

-22.225 

2.485,  22.233 

8 

.177,  26.915 

2.301, 

-25.346 

2.508,  25.353 

9 

.159,  30.103 

2.349, 

-23.470 

2.533,  28.477 

10 

.144,  33.283 

2.393, 

-31.598 

2.558,  31.603 

11 

.132,  36.458 

2.433, 

-34.727 

2.583,  34.732 

12 

.121,  39.628 

2.469, 

-37.859 

2.606,  37.863 

13 

.112,  42.793 

2.501, 

-40.992 

2.628,  40.996 

14 

.104,  45.955 

2.529, 

-44.127 

2.648,  44.130 

15 

.096,  49.114 

2.554, 

-47.263 

2.665,  47.266 

16 

.090,  52.270 

2.575, 

-50.401 

2.680,  50.403 

17 

.084,  55.424 

2.593, 

-53.539 

2.693,  53.541 

18 

.078,  58.576 

2.608, 

-56.679 

2.703,  56.680 

19 

.073,  61.725 

2.619, 

-59.820 

2.710,  59.821 

20 

.069,  64.872 

2.627, 

-62.962 

2.714,  62.962 

21 

.064,  68.017 

2.631, 

-66.105 

2.715,  66.105 

22 

.060,71.160 

2.630, 

-69.250 

2.712,  69.249 

23 

.055,  74.300 

2.626, 

-72.396 

2.705,  72.394 

24 

.051,  77.437 

2.616, 

-75.544 

2.694,  75.542 

|  Table  3.  Auxiliary  Functions  for  the  KGSB  and  HSGB  Solutions 

Wr)  =  Z\k  +  difZji ,  +  e*Zi3*  4-  /i-Zig* 

Z?kB(r)  =  Zik  +  dkZgk  +  CkZuk  +  fk  Z'iQk 
ZzkB(r)  —  Zzk  +  dkZgk  +  CkZ\Sk  +  f*Zi\k 
Z*kB (r)  =  Z<*  +  dkZxoi,  +  e*Zi6*  +  fkZnk 
ZskB(r)  —  Zsk  +  dkZnk  +  tkZnk  +  fkZizk 

ZekB{r)  =  Zek  +  dkZ  12*  -f  tkZ\ik  +  fkZuk 

The  three  sets  (Z,*),* r...n,  (Zjk)jm  13...1S*  (Z,*);*i«...24 
satisfy  the  homogeneous  form  of  the  boundary  conditions 
satisfied  by  the  set  (Zjfc)j=i...e.  At  the  interface: 

for  j  =  7  13  19 

Z]k(ri)  =  1  0  0 

Zo+i)fc(ri)=0  1  0 

Z(:+z)k(r\)  =  0  0  1 
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Table  4.  Liquid  Moment  Coefficients 


Cmp/  =  *(1  -  «)  {[/  -  2  -  co Z«,(1)M73  -  (l/2A)Ec*P*Ztt(l)} 

Cmp<  =  i(l  -  ie)  [(2  -  f)Bi  +  Co/r*  qlrZeo(r)r3dr  +  (1  /A)Zckgkf?0  ?irZ6i(r)r3dr] 

Cmpr  =  -J9ir0(l  -  i<)  {[(/ -  2)r0  -  c0Zeo(r0))A3 /Z  -  (1/24)  £  CfcDfcZ64(ro)} 

Cmui  =  ^  {co[iZ»(l)  -  (A2/3)Z3o(1)]  +  Ec* 

Cmu<  =  {-2/3,  +  co/,1,  93r[2Zi0(r)  -  ZM(r)]r  dr  +  T.ck\kh„^  q3r\2Zlk(r)  -  Z2*(r)]r  dr} 

Cmvr  =  {c0[ir0Zso(r0)  -  (A3/3)Z3o(r0)]  +  Eek 

where  all  summations  are  for  k  =  1  to  N  and  where 
Bi  =  [1  -  (1  -<?i)rj  -  ?ir£]/4 
B2  =  [1  -  (1  -  q3)r\  -  93ro]/2 
c0  =  /  for  HSG  or  HSGB 
=  0  otherwise 
Di,  —  f*A  xGk(x)dx 
qi r  -  qi  (r  <  r,) 

=  1  (r  >  r,) 

Kr  =  <73  =  jj?  (r  <  r,) 

=  Mr  >  r,) 

7i  =  <73  =  1  for  a  single  liquid 
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Table  5.  Boundary  Layer  Modifications  to  tiie  Liquid  Moment  Coefficients 


For  SW,  SWB,  KGS  and  IvGSB: 

Cmvt  =  ^  {-(1  -  3A)2fB,  +  ,  ck(Xkhk  -  3gk)fr\  q3r[2Zlk(r )  -  Zik{r)\r  dr} 

For  SW  and  SWB: 

In  Cmpi,  replace  Z6*(  1)  by  Ze*(l)  -  2i6ae7Fl(\) 

In  CmpT,  replace  Z6*(r0)  by  Z6k(r0)  +  2iSt,e2Fl(r0) 

In  Cmvi,  replace  Z5*(  1)  by  Zik(  1)  -  $~le2Z4*(l) 
and  replace  Z3k(  1)  by  Z3k(l)  -  i6~ie1Fi(  1) 

In  CmvT ,  replace  Z6*(r0)  by  Z5i(r0)  +  <5^ e2Z4t(r0) 

?.nd  replace  Z3*(r0)  by  ;3*(fo)  + 

where  3  ~  (1  +  *)\/0  +  /)/?e„/2 
F\(r)  =  Z2k(r)  —  Z,*(r)  +  r 

h  =  M?" 17  2 
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Table  6.  Components  of  the  Liquid  Angular  Momentum 


L  =  L\exc  +  Lieyc  +  L^iic 
where 

Lj  =  a4  fo1  fl4A  fo*  p(f)N,r  dj>  dx  df 
and  where 

p(f)  =  0,  0  <f<f/(x,0) 

=  Pa,  r/(*,  0)  <  f  <  f,(x, t£) 

=  Pt,  f,(x,  0)  <  f  <  1 
iV,  =  -yrVg  4-  K{(xVr  -  rV'x)sin^>  4-  xVj  cos^>] 
jV2  =  KrVg  —  7[(xVr  —  r!4)sint/>  4-  xV*  cos  V’] 

A'3  =  (xV”r  —  r!4)  cos  ^  —  xVg  sin  0 
For  small  K, 

Lx  =  Imia^^Bx 
Li  =  m£,a2dA'f?{iS} 

X3  =  m£,a2<£A'f?{S} 
where 

5  =  /  +  (1  —  ?i)/Fj(r1)  +  ?i?4/Fj(ro) 

/  =  —(1/2 A)  Jr‘0  gir[f(»r  -  in)  +  x(jt  +  iii  4-  2tx]f  rfrdx 

/Fj(r)  =  +  *(1  —  /)x]x  dx 

^4  =  1  for  partial  fill;  zero  otherwise 


Table  7.  Liquid  Moment  Coefficients  Computed  from  the  Angular  Momentum 
For  HSG,  HSGB,  S  and  SB: 

Cmm  =  (1  -  ie){2iDl  +  /[F3  +  (1  -  9i)/F2(ri)  +  <?i (/  +  <74  -  l)fa(r0)]} 

For  KGS  and  KGSB,  add: 

(ACmm)j  =  -(1  -  «)/  {[(A/3  -  1)//32]F4  -  (1  -  /)F2(1)  -  (1  -  <7i)F2(r,)  -  qiq*F2(r0)\ 
For  SW,  add  (ACfnm)l  plus 
(ACmm)j  =  —(1  -  ie)f6aEFb(\)e2 

For  a  central  rod,  using  SW,  add  (A Cmm)v  (ACmm)2  and 
(ACmm)3  =  -(1  -  ie)f6bEFs(r0)e7 
All  summations  are  for  k  =  1  to  N  and 
Bz  =  [1  -  (1  -  <?,)r2  -  qxrl\/ 2 

E  =  1  -  e(r°-I)/*. 

F2( r)  =  ^  {(2A2/3)  [/  -  CoZ,0(r)]  -  (1/A)  £  c*DkZ,*(r)} 

Fa(r)  -  -5jf=7J  [2/  -  (UA)Zck9k[2Zu(r)  -  Z»(r)j]  -  ^  £ctff*Z»(r)} 

F3  =  -»(/?!  +  2£?3A2/3)  -  co/0  +  E  c*A(A*  -  2gk/A\k) 

Fi  =  2fB,-(\/A)Zc3gJk 

Fs(r)  =  r  {/.42/3  +  (1/2 A)  £  c*{Z?*[Zj*(r)  -  Zi*(r)]  —  2r<7*Z.u(r)/A* } } 

2*  =  Jr‘#  <7uZ4*(r)r2  dr 
U  -  Ir\  g«r[2Z,*(r)  -  Z2k(r)]r  dr 

a  =  -  /)/( 1  +  /) 

<  ,  -1/3 
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Table  8.  Liquid  Moment  Coefficient  Differences  for  Partial  Fill  and  L  0 

A OnW/L2  =  -i(l  -  «<)[/  -  Z«,(l)  +  (1/it)  2:  A*  »e«Zi*(l) 

ACmW/I5  =  (1-  «)(/2e)-‘[Zu(l)  -  (1M)E  A^cu^l) 

A Cmpr/L7  =  »r0(l  —  *«)(/ro  “  2w(*“o)  +  (1/-4)E  AJ*Ci*Z8*(ro) 
ACmw/Ls  *  — r0(  1  -  M)(*e)-'(Zu(ro)  -  (1M)E  A?  Ws*(ro) 
where  all  summations  are  for  k  »  1  to  N 
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Figure  2.  (by  HSG,  volume  and  surface  method 

(Re  *  500.  A  »  3.  0.1). 
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Her  data. 
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Appendix  A:  Effect  of  Cenfcer-of-M&ss  Axial  Offset 

If  the  projectile’s  center  of  mass  is  not  at  the  center  of  the  cylinder,  the  liquid  moment 
should  be  computed  relative  to  the  projectile’s  center  of  mass.  Since  the  projectile’s  coning 
motion  is  relative  to  its  center  of  mass,  the  boundary  conditions  on  the  liquid  are  also 
changed.  In  this  appendix,  we  will  consider  both  of  these  effects. 

We  will  locate  the  projectile’s  center  of  mass  relative  to  the  center  of  the  cylinder  by 
the  dimensionless  distance  L  —  h/a.  This  distance  will  be  negative  if  the  projectile  center 
of  mass  is  to  the  rear  of  the  cylinder  center.  Eqs.(10-12)  become 


x  —  L  =  x  -  L  +  fR{Kt-^}  (41) 

r  =  f  -  (x  -  L)R{Ke~i*}  (42) 

0-0  =  (43) 

The  new  boundary  conditions  become: 

(a)  on  the  cylindrical  wall  (f  =  1): 

=  *(1  -/)  (44) 

li(l,*)  =  — *(  1  —  f)(x  —  L)  (45) 

!i(l,x)  =  -(1 -/)(*-£)  (46) 

(b)  on  the  end-walls  ( x  ±  4): 

Ul(r,±4)  =  i(  1  —  f)r  (47) 

tt(r,  ±4)  =  i(l-/)(T4-I)  (48) 

t(r,±4)  =  (1-/)(TA-X)  (40) 


If  we  add  a  subscript  ‘c’  to  and  Cp  to  denote  those  variables  for  L  =  0,  then 

more  general  variables  can  be  defii  :i  by  the  relations: 


UL 

1L 

li  ~  W 

a 

c. 


uic  +  W 
tie  - \  L 


N 


xZil(r)  +  '£icLi,Z<k(r)Gk'T) 


km\ 


iV 


«(1  -  /)  +  fZXL(r)  +  /  £  rLkZlk(r)Hk(x) 


km  l 


lie  ""  *lie  4’  Lf 


N 


ZtL(r)  +  Y,c"Zn(r)Hk(x) 


km  l 


2e  +  L 


H 


(1  -  f1)r  4-  /Z«t(r)  +  /  £  cuZM(r)ffk(x) 


kml 


Cp.  +  Lf 


s 


-fr  +  Z.i(r)  +  £  ctkZ6k(r)ffk(x) 
*■» 


(410) 

(411) 

(412) 

(413) 

(414) 


S3 


The  terms  involving  Zj£,  must  satisfy  the  linear  Navier-Stokes  equations.  This  means  that 
they  are  solutions  of  Eqs. (52-57)  when  index  k  is  replaced  by  the  symbol  L  in  each  subscript 
and  when 

\i  =  0,  At!  =  1,  \l2  —  0  (-415) 

It  can  easily  be  seen  that  the  other  terms  multiplied  by  L  and  outside  the  summation  are 
also  solutions  of  Eqs. (52-57)  with  k  replaced  by  L. 

We  now  require  that  ti.  :  terms  multiplied  by  L  satisfy  the  boundary  conditions  mul¬ 
tiplied  by  L.  For  r  =  1, 

Z1L(1)  =  Z2L(l)  =  Z4L(1)  =  0  (A16) 

For  the  interior  boundary  condition0: 

(a)  for  a  fully-filled  cylinder: 

ZeU  0)  =  Z2L(  0)  =  Z4l(  0)  =  o  (A17) 

(b)  for  a  partially-fi'led  cylinder: 

o.  ,  s  .  *roZ,t(ro)  2 Z[L(r0)  , 

Z«(r,)+  — 7 - =  / 

ZiJr0)  =  0 

r0Z3L(ro)  ~  iZ2L(ro)  =  0 

(c)  for  a  cylinder  with  a  central  rod: 

Zu(ro)  =  Z2L(ro)  -  Z4L(r0)  -  0  (-421) 

Note  that  Eqs. (55-56)  contain  only  Z4l  and  Z%i.  Since  the  boundary  conditions  for 
these  variables  are  homogeneous, 

Z4L(r)  =  ZiL(r)  =  0  (.422) 

On  the  end-walls: 

.v 

£c,,*Z«(r)  =  0 
,v 

Z\dr)  +  ^cLkhkZu,(r)  =  0 

km\ 

V 

Zu.(r)  +  £  c,.kUkZ2k(r)  =  0 

ka  [ 

lor  the  fully  filled  and  central  rod  cases,  all  the  boundary  conditions  are  homogeneous 
and  therefoie 

Z,ar)a  0  (.420) 

M 


(.423) 

(.424) 

(.125) 


(-41S) 

(-419) 

(A20) 


For  these  cases,  Eqs.(A23-A25)  can  be  satisfied  by 

cLk  =  0  (-427) 

and  only  the  pressure  components  of  the  liquid  moment  are  affected.  If  moments  are 
computed  with  respect  to  the  actual  center  of  mass,  the  entries  for  Cmvi  and  Cmpr  in  Table 
4  must  be  modified  by  the  following  increments: 

A Cmpl  =  -*(1  -  ie)fL2  (428) 

A Cmvr  =  i>o(l  -  ie)JL 2  (A29) 

For  the  partially-filled  case,  however,  the  presence  of  the  free  surface  makes  the  situa¬ 
tion  much  more  complicated.  Four  of  the  entries  in  Table  4  must  be  modified;  the  required 
increments  are  given  in  Table  8. 
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Appendix  B:  Linear  Moment  from  Angular  Momentum 


If  we  apply  Eqs. (10-12, 19-21)  to  the  definitions  of  the  Nj's  (Table  6),  we  obtain  the 
following  linear  expansions: 

Ni  =  a4>[f2  +  Kgi]  ( Bl ) 

JV2  =  -a<j>  rx cos rj> +  (K/2)R{iN  +  g2}]  ( B2 ) 

iV3  =  -a<t>  rising  +  (K/2)R{N  +  53}]  ( BZ ) 

where 

N  =  r(ir  —  m) -r  x(u  +  in  f- 2t'x) 

9j  = 
ri* 

j  9jdx})  =  0 

Since 


p(r)  =  0,  0  <r<f/(i,0) 
=  Pit  r/ixtfy  <  r  < 
=  /h,  fj(i,0)  <  r  <  1 


we  can  write 


1  f1  fA  f2* 

P\  JQ  J_J0  PNS  ^  dxdf  =  MX}  -f  +  qxMzi  +  M4j 


where 


r\  fA  rl* 

»  /  /  /  quNjf  dxf)  dx  df 

<f»  J—A  JO 

A/?,,  M3],  MA}  =*  J  a  Jq  J  Njrdfdxfrdx 

(rAtrB)  m  (f/,f0)  for  M7j 

*  for  M3i 

■  for  MAi 

For  a  central  rod,  f/  =*  f0  and  A/j,  =»  0.  Linear  expansions  of  the  Af,-/s  yield 

Lx  m  2mLa24>Bx  (£5) 

Lq  *  mtaa^A'if{«5}  (S6) 

Ls  m  mLa2j>KR{S }  (£7) 

where 


S  =  I  +  (1  -  q\)fF3(r\)  +  qxqifFiiro) 

I  *  -(1/2A)  J  j  ^qirNrdxdr 

fFi(r)  m  -  4(r,i)  +  i(l  -  f)x]x  dx 
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and  where  =  1  for  the  partially  filled  case  and  zero  otherwise. 

Note  from  the  definition  of  .V  that  the  Mt] 's  involve  the  volume  integral  Vi  of  f(ir  —  m) 
and  the  volume  integral  V2  of  i(y.  +  «li  +  2ix).  We  now  derive  a  relation  between  V2  and  V2 

in  order  to  eliminate  V2.  The  continuity  equation  for  the  perturbation  velocity  functions 

isi'io,n]; 

d(fn)  _  Dw 

C  =  — r- - ill  4-  r  —  =  0  ( B8 ) 

or  ar 

If  we  integrate  by  parts  the  volume  integral  of  qXrfxC.  we  obtain 


La{  M!,  +  i  di  +  [aqir[iui  }1.Ar2  df 

~  /  +  »ll)i  +  mf}r  dx  dr  =  0 

*  *  —  A 

Adding  ix  and  if  terms  as  needed  to  obtain  the  desired  expressions,  we  have 

J  A  |  +  ,i)]M  +  [?i^2(u  +  l-Mro}  *  +  J  9ir  f2  df 

—  I  /  9ir((tt  +  «1£  +  2ix)x  +  (u>  —  ir)r]r  dx  df  =  0 
Jr0  J-A 

which  reduces  to  the  desired  relation: 

(1/2 A)  J  y  ^  9ir(li  +  Jil  +  2ix)xf  dx  dr  =  (1/2.4)  ^  J  q,T(if  —  uf)f2  dx  df 

+0  -  f)fhiF2(r0)  +  (1  -  <]i)Fi(ri)  —  F2(  1)] 
+  «/B2(2.42/3)  -  i/Z?,  -  fG 

where 


(BIO) 


(511) 


2 

Fj(r)  =  ofl”-"/)  {("‘^2/^)  “  co^«°(r)]  “  (1/-4)  c*5*Zj*(r)} 

fG  =  f  q\,[i(l  -  f)r  -  UL(fsA)]f7  dr 

Jr0 

According  to  the  boundary  conditions,  Fj ( 1)  and  G  are  zero.  This,  of  course,  may  not  be 
the  case  for  a  finite  number  of  eigenfunctions  or  for  boundary  layer  solutions  based  on  an 
approximate  form  of  the  continuity  equation.  However,  F}{  1 )  =  0  is  a  valid  assumption 
for  HSG  and  HSGB.  The  actual  nonzero  value  of  Fj(  1)  will  be  computed  for  the  boundary 
layer  theories.  G  will  be  computed  exactly  in  all  the  theories. 

F.q  (Bll )  ran  now  be  used  to  derive  the  formulas  in  Table  7. 


Appendix  C:  End- Wall  Boundary  Layer  Relations 


The  usual  boundary  layer  assumptions  can  be  made  for  the  linearized  Navier-Stokes 
equations  [Eqs. (29-32)]  in  the  vicinity  of  the  end  walls  (x  =  ±,4).  The  subscript  “be”  will 
denote  end-wall  boundary  layer  functions. 


dine  ,  dm 

-^7  +  i^-^  +  r— 


=  0 


*'(/  -  l)lZte  +  2 iihe  = 

^  = 

dx 


1  d2^ 
Re  dx 2 

Re  dx 3 

0 


(Cl) 

(C2) 

(C3) 

(C4) 


If  the  longitudinal  symmetry  of  the  boundary  conditions  is  used,  the  solutions  and  u*, 
to  these  equations  take  simple  forms  closely  related  to  those  given  in  References  [4]  and 
[12]. 


It*  =  (,7/2)[-^(r)C.(r)  +  B(r)C4(x)] 
Ife  =  {f/2)[A{r)Ga{x)  +  B{r)G,(x)\ 


(C5) 

(C6) 


where 


G„  =  sin(tax)/  sin(ia>l) 

C\  —  sin  (i/dx)/  sin(«/M) 

a  =  a  -  i)((3  -  /)«e/2J,‘,a 
fi  =  (1 +>)[(! +/)B«/2]‘/J 

The  functions  di(r)  and  B(r)  are  determined  by  the  no-slip  conditions  at  x 
denote  the  left  sides  of  Eqs. (33-36)  by  w,  u,  v,p,  we  can  write 

Si(r,  A)  =  (u  +  =  —  *(1  -  f)A 

u(r,/l)  a*  (v  +  =*  “(1  -  f)A 


Thus 


f  Mr)  -  -(»  +  iu)„4  at  -»/£c*Zj*(r) 

/B(r)  =  -{v  -  tu  —  2(1  —  /)x]*.4 

=  -«/{2/\4-£c*(2Zlt(r)-Z,*(r)l} 

The  continuity  equation  can  now  be  used  to  obtain  a  relation  for  w^: 


Uiu 


-1 


20-/) 


(1  +/)C.(z)  (3-/)C»(x) 

r  ^ 


5ui(rM) 


5x 


The  boundary  condition  for  the  velocity  normal  to  the  end  walls  is 
ui(r,  ±A)  a  ( u>  +  Wu  )..±a  »  *(  1  -  f)r 


A.  If  we 

(C7) 

(C8) 

(C9) 

(CIO) 

(Cll) 

(012) 


or 


Y^ck24k(r)[cos(.\kA)  +  A k6c  s:n(At.4)]  =  0 


(C 13) 


where 


<5,  = 


-1 


2(1-/) 


1  +  /  .  3-/ 

a  t 


-(!+«') 


2(1  -fW2Rc 

Thus  Ak  must  be  selected  so  that 

cos(AkJ4)  +  Xk6e  sin(Ak,4)  =  0 


1  +  /  +«(3-/) 


\/3  -  /  n/TT7 


(C 14) 


The  boundary  condition  for  the  normal  velocity  at  the  lateral  wall  is  affected  by  the 
presence  of  the  boundary  layers  at  the  end  walls. 

li(l,x)  =  u(l,x)  +  m,(l,x)  =  -«(1  -  f)x  (05) 

This  reduces  to 

^ckZlk(l)Gt(x)  +  (.72)(-^(l)G.(x)  +  i?(l)Ck(x})  =  /x  (05) 

Now 

^(1)  »  0  (07) 

i?(l)  =  -2i/(A  -  ]Ta*$*)  (08) 

where 

9k  =*  Gk(A) 

We  select  the  Zu’s  to  be  unity  and  the  ak’s  to  be  the  least  squares  fit  of 

H  ak[Ck(x)  -  ?kC»(x)]  =■  x  -  AG>(x)  (09) 

Eq.(ClO)  differs  slightly  from  the  usual  conditions  for  the  aks.  For  very  low  Reynolds 
numbers,  this  difference  is  significant  and  hence  (09)  is  used  in  the  KGS  and  SW  codes. 
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Appendix  D:  Lateral- Wall  Boundary  Layer  Relations 


From  Eqs.(29-32),  the  boundary  layer  equations  near  the  lateral  wall  (r  =  1)  can 
easily  be  written.  The  subscript  ubF  will  denote  lateral  wall  boundary  layer  functions. 


dr 

dJ*L 

dr 

*(/  -  Ikw 

*(/  -  l)l£w 


1 

Re  dr> 
i 

Ye  dr  3 


The  solution  of  Eqs.(D3-D4)  can  be  written  in  the  form 

Uu  »  u,(x)erp{(r  -  l)/£«)  +  t*(x)exp{(ro  -  r)/S. } 
UU4  “  «i(x)exp{(r  -  l)/£.)  +  w<>(x)exp{(r®  -  r)/<$„] 

The  no-slip  condition  at  r  sc  1  yields 


(D 1) 
(D2) 
(D3) 
(D  4) 

(D5) 

(D6) 


*  t(l,x)  +  Vj  +«|tt8  *  -(1  -  /)x  (D7) 

it(l,x)  »  iir(l» x)  +  u?t  +  ei<i>Q  ■  t(l  -  /)  ( £>8 ) 

where 


«i  =»  <xp((r0  -  1)/*.) 

1£  •  tp  +  Ifi*, 

For  the  partial-fill  case,  no  shear  at  the  free  surface  yields 


Thu. 


«i*h  -  u.  ■  0 

f,W,  -  10Q  m  0 


zilzilfLlillifl. 

1+ef 

i(l  -  /)  -  ft(l, x) 

1  +  e?  1 


Uo  *  «|Vj 


U’o  m  e|tl>} 


For  a  central  rod, 

£(r0,x)  +  «iVj  +  Vo  m 
lii(pOiX)  +  eiu>|  +  wq  m  *(  1  —  /)r0 


(DO) 

(DIO) 

(DU) 

(DU) 


(£>13) 


Thus 


V\ 

I'o 

Wj 

«’o 


~(1  -  /)*  -  i(U)  _  ei[r(l,  x)  -  r(r0,x)] 

1  +  t\  1  -  ej 

-( 1  -  f)x  -  r(rQ,x)  £i  [£:(!,  x)  -  r(r0,  x)] 

1  +  e,  1  -  ej 

i(l  -  /)  -  1£(1,  x)  _  ei[li:(l.  x)  -  lh(r0,x)] 

1  -f  ej  1  —  ef 

i(l  -  f)rQ  -  tr(r0.  x)  et[ir(l.  x)  -  »l:(r0,  x)] 
1  +  e,  1  —  e. 


The  continuity  equation  can  now  be  used  to  obtain  a  relation  for 
ni*(r,x)  =  <5a{u,exp{(r  -  l)/6a]  -  u0c^p[(r0  -  r)/<$a]} 

where 

6,  =  (1  +  »)[2(1  -f)Re}-'/7 

U]  =s  U),  —  U’j 
uo  —  it’  o  —  “ott'o 

For  a  partial  fill, 


ui  = 


1  +  «? 


u0  =  e,u, 


where 


For  a  central  rod. 


■•••hr 
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<Xrjl) 
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1  +  r, 
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'<«  e,A* 
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Thus 
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(D17) 

(D1S) 

(D19) 

(D20) 

(021) 

(D22) 

(023) 


(024) 


(025) 


02 


uM(l,r)  = 

partial  fill 

(D26a) 

= 

rod 

(D26b) 

uu(r0,x)  = 

0 

partial  fill 

(D27a) 

= 

(i+ Ju‘+( 

ft)* 

rod 

(D27b) 

The  A*  terras  in  the  rod  expressions  above  are  usually  much  smaller  in  magnitude  than 
the  other  terms  and  will  be  neglected. 

The  boundary  condition  at  r  =  1  can  be  written  as 

«(l,x)  +  «*e(l,jr)  +  u*(l,r)  =  —  i(l  -  f)x  (Z>28) 


or 


N  (  <5ae2  \  du(r,x 


=  ~«(1  ~  /)*  - 


r*l 


where 


1  —  e* 

n  =  1  for  a  central  rod ;  2  otherwise 

1  +  e" 


The  boundary  condition  at  the  rod  surface  (r  =  ro)  is  similar  to  (D29): 


u(r0,l)  + 


(  Str0e j  \ 

0u(r,x)' 

V  ro  +  / 

dr 

=  — t(l  -  f)x  —  Uu(r0,x) 


(D29) 


{D  30) 
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Intentionally  left  blank. 
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a 


List  of  Symbols 

radius  of  the  liquid-payload  cylinder 


ai  coefficients  in  a  least-squares  fit  of  x  to  a  sine  series 

(a*  =  0  for  k  =  0,2,4,—) 

ak,  bk  coefficients  in  the  definitions  of  Z\kB  and  Eqs.(122-123) 

A  c/a,  fineness  (aspect)  ratio  of  the  cylinder 

b  radius  of  the  inner  cylindrical  free  surface  in  a 

partially-filled  cylinder 

3  r*1  +  AJ  -  i(l  -  f)Re  [replace  Re  by  for  the  two-liquid  case] 

Bx  [l-(l-$,)rf-?ir*]/4 

Bt  [1  -  (1  -  $3>?  -  <hrl\n 

B3  [l-(l-fll)rf-«»r|]/2 

c  half-height  of  the  liquio-payload  cylinder 

c<,  /  for  HSG  or  HSGB;  »  0  otherwise 

c*  .onstant  multipliers  in  the  definitions  of  <?„*,  Eqs.(46-51);  Fourier 

coefficients  choeen  to  satisfy  boundary  conditions 

cmi,„cmpi,cmpr  that  part  erf  Cistd  +  iCuu  due  to  thj  yressere  on  the  end-walls, 
lateral  wall,  and  rod,  respectively 

cmwT  that  part  of  Ccsu  4  iCu»4  due  to  the  vitcont  strese  on  the  end-walls, 
lateral  wall,  and  rod,  respectively 

CFD  Computational  Fluid  Dynamics 

Ccsu  +  iCuM  (A/*  +  iMn)/mia*i3TK,  the  liquid  side  moment  and 
in-plane  moment  coefficients 

Cunu  Af#</"*tV 'i*rK3,  the  liquid  roll  moment  coefficient 


C9  pressuie  coefficient,  Eq. (150) 

d  radius  of  a  central  rod  within  the  cylinder 


Dk 

t 

ex,  ev, ej 

CXC  '  '  £  TC 

e\ 

«2 

«3 

e< 

c$ 

c« 

f 

i 

h 

Ok 

Gk 

h 

hk 

lh 

HSG 

hsgb 


fl\  xGk(x)dx 

(R/ Rq)1/2 ,  error  measure  of  the  fit  to  Eqs.(101-103) 
unit  vectors  in  the  earth-fixed  cylindrical  system 
unit  vectors  in  the  earth-fixed  rectangular  system 
unit  vectors  in  the  coning  aeroballistic  system 
exp[(r0  -  l)/<$a] 

(1  —  e")/(l  +  e")  where  n  =  1  for  a 
central  rod;  otherwise,  n  =  2 

eip[(r,  -  l)/<5„] 

«72/J[l  +  +  (1  "  oM]  -  (1  -<?,)(!  +  el)K 

[1  +  <75  -  (1  -  +  <?5  +  (1  -  <75)^3] 


[1  +  ?S  +  (1  -  1  ~ 

(1  -«)r 

_2j(l -/)/(!+/) 

[(3  -  /)(!  +  /)]*^/(l  -  /) 

Gk(A),  see  Eqs.(  101-102) 

function*  of  x  in  the  expres*ion*  for  u,  i>,  and  g,  Eq*.(34-3C);  Get 

distance  from  the  projectile  center  of  mas* 
to  the  cylinder  center  of  mass  (positive  if  the 
projectile  c.m.  if  forward  of  the  cylinder  c.m.) 

Hl(  A ),  see  Eq.(  103) 

function*  of  1  in  the  expression  for  in,  Eq.(33);  /f0  =  1 
Hall-Sedney-Gerber  method.  Ref.  [2 2 J 
HSG  Binary  method 


Ik 


Jr\  quZ<k(r)r2  dr 


Jo,Jx 


k 

K 

KGS 

KGSB 

A'* 

A'o 

L 

L 


mi 

M 

Af*«,  A/yc ,  A/Jc 

P 

Po 

p(r,r} 
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7< 

7$ 


complex  Bessel  functions  of  the  first  kind,  order  0  and  1 
(with  argument  /^Atr) 

summation  index  for  the  perturbation  variables, 

Eqs.(33-36),  k  =  0,1,2,...,  N 

sin  a, 

Kitchens-Gerber-Sedney  method,  Ref.  [10] 

KGS  Binary  method 

Ke~'*  =  Koc'u*~t) 

the  value  of  A'  at  it<$>  =  0;  A'  =  A'0e‘r * 
h/a 

angular  momentum  of  the  liquid 
2  xa2cpL 
liquid  moment 

components  of  the  liquid  moment  in  the  coning  system 

liquid  pressure 

value  of  p  at  r  =  0,  A'  =*  0 

complex  pressure  perturbation,  nondimensionalized  by  pia 3<j>3 
Pi/P\ 

i/  >  /  t/j 

/i,/p,  =  7,/7j 

1  for  the  partial  fill  case;  0  otherwise 

7./^’ 
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q\T  ?i  (for  r0  <  r  <  or  1  (for  rj  <  r) 

q2r  q2  (for  r0  <  r  <  rj )  or  1  (for  rx  <  r) 

93r  93  (for  r0  <  r  <  ri)  or  1  (for  r2  <  r) 

r,  r  earth-fixed  and  coning  radial  coordinates,  respectively, 

nondimensionalized  by  the  radius  a 

rc  earth-fixed  radial  coordinate  of  the  cylindrical  wall, 

nondimensionalized  by  a 

Tf,ff  earth- fixed  and  coning  radial  coordinates  of  the  free  surface  in  a 

partially-filled  cylinder,  nondimensionalized  by  a 

ri,fi  earth- fixed  and  coning  radial  coordinates  of  the  two-liquid  interface 

(for  small- amplitude  coning  motion),  nondimensionalized  by  a 

rT,fr  earth- fixed  and  coning  radial  coordinates  of  the 

rod,  nondimensionalized  by  a 

r0  (a)  for  a  fully-filled  cylinder,  0 

(b)  for  a  partially- filled  cylinder,  the  radius  b/a  of  the  inner 
cylindrical  surface  (for  no  coning  motion) 

(c)  for  a  cylinder  with  a  central  rod,  the  rod  radius 
nondimensionalized  by  a 

r»  earth-fixed  radial  coordinate  of  the  two-liquid  interface  (for 

no  coning  motion),  nondimensionalized  by  a 

(rf )  inner-liquid  value  at  r1 

( r *)  outer-liquid  value  at  rj 

R  function  minimized  in  the  determination  of  the  c*’s,  Eq.(104) 

R  position  vector 

R{  }  real  part  of  a  complex  quantity 

Re  Reynolds  number,  (or,  for  two  liquids, 

the  outer-liquid  Reynolds  number) 

Re%  Re  for  v  m  v% 


M 


/ 


R o 

Rl,R-2,  R3 

s 

SB 

SW 

SWB 


value  of  R  for  c*  =  0 

functions  used  in  determining  the  c^’s,  Eqs. (101-103) 

Stewartson  inviscid  method 

Stewart  son  Binary  inviscid  method 

Stewartson- Wedemeyer  high- Re  method 

Stewartson- Wedemeyer  Binary  high- Re  method 

complex  radial  velocity  perturbation, 
nondimensionalized  by  a4> 


v(r,x) 


V 

Vt,Vr,Vt 

m(r,x) 


x,x 


complex  azimuthal  velocity  perturbation, 
nondimensionalized  by  a<fr 

velocity  vector 

velocity  components  in  the  earth-fixed  cylindrical  system 

complex  axial  velocity  perturbation, 
nondimensionalized  by  a<fr 

earth-fixed  and  coning  axial  coordinates,  respectively, 
nondimensionalized  by  the  radius  s 


X  axis  along  the  projectile’s  axis  of  symmetry 

X «,  Y„  Z,  axes  in  the  earth-fixed  system,  X,  initially  along 

the  velocity  vector,  Z,  downward 


r0,r, 


complex  Bessel  functions  of  the  second  kind,  order  0  and  1 
(with  argument  mA*r 


Znk 


complex  perturbation  variables  formed  from  (t&*,ti*,  t>*), 

Eqs.(46-51) 


a  (1  -  0^(3  -  /)/&/ 2  -  -1^(3  -/)/(!+/) 

at  total  angle  of  attack;  the  angle  between  the 

X  and  X,  axes 


99 


<3  (l  +  i)y/{l  +  f)Re/2 


7 

<5co 


AZ„*b 


A* 

At  i 
Air2 
Ml-/1! 

1/ 


l/i 


r/-* 


/v. 


■  />? 

r 


coso( 

1-H 

r  —  1/2 

M2 

-(>+<)  W  ,  <(3-/) 

2<i-/)vTfli  1^5-7  t  yrr? 

<5C  for  iZe  =  Ac* 

ZnkB(rt)  ~  ZnkBirD,  the  jump  in  Z„*b  at  the  two-liquid  interface 
(yaw  damping  rate)/(coning  rate)  =  (A* /K)/(<j>a) 
earth-Exed  azimuthal  coordinate 

complex  eigenfrequencies  arising  in  the  expressions  for  the  perturbation 
variables,  Eqs. (33-36),  where  A0  =  0 

0  (if  k  =  0),  A*  (if  k  >  0) 

1  (ifk  =  0),  A*  (if  k  >  0) 

outer-  and  inner-liquid  dynamic  viscosities,  respectively  (p}  =*  u}p}) 
kinematic  viscosity  of  the  liquid 

average  two-liquid  value  of  the  kinematic  viscosity,  Eq.(126) 

outer-  and  inner-liquid  kinematic  viscosities,  respectively 

density  of  the  liquid  (or,  for  two  liquids,  the 
density  p\  of  the  outer  liquid) 

outer-  and  inner-liquid  densities,  respectively 

<?t ,  the  roll  angle 

spin  rate  (assumed  positive  and  constant);  the  axial  component 
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of  the  projectile’s  angular  velocity  with  respect  to  an  inertial  frame 

<j>a  the  angle  between  the  Ze-axis  and  the 

normal  to  the  X  —  Xt  plane 

<j>a  precessional  rate  of  the  angular  motion  performed 

by  the  projectile  about  the  velocity  vector 

rp  9  —  <pa  =  9  —  r<p 

tp  coning  system  azimuthal  coordinate 

P  angular  velocity  of  the  coning  aeroballistic  system 

Misc.: 

(')  d(  )/dt 

(  Y  d{  )/dr 
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This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments/answers  to  the  iteas/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Suaber _ TR-3074 _ _Date  of  Report  ^ 
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etc?  If  so,  please  elaborate. 
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